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ABSTRACT 



The theory of singular lines on the Parameter Plane 
and in Parameter Space are derived, and applied to control 
system design. 

A method for design of compensation of linear feed- 
back control systems using singular line theory is 
presented, and application of this design method for 
self-adaptive control systems is considered. 

General design steps and required procedures for using 
this design method are summarized and examples are pres- 
ented. 
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I. INTRODUCTION 



Prior to 1959 analysis and design of control systems 
was primarily performed by classical methods. These 
methods are adequate when there is a single variable 
parameter and a single loop system. 

In modern control engineering, where multiloop systems 
with characteristic equations with two or more variable 
parameters usually result, the above mentioned design 
methods become in general very complex. The designer has 
to complete several designs in a trial and error process 
and to choose the most suitable solution and to make 
usually some arbitrary choices based on his experience. 

Methods for studying the effect of variations in 
parameters of a multivariable, multiloop control system 
have been developed by Mitrovic [1] , Siljak [2] , Thaler 

[3] , and other investigators. All these methods are based 
on the same concept and their dependence on a computer is 
characteristic . 

Mitrovic' s [l] method which was introduced in 1958, 
allowed study of the effect of variations of two variable 
parameters on the location of the characteristic roots of 
a system. This method consists of specifying as variables 
the two lowest order coefficients of the characteristic 
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equation. Using then the characteristic equation as a 
mapping function, constant zeta, omega and sigma curves 
are transformed from the s-plane into the variable co- 
efficient plane. This method permits the designer to 
adjust these variable coefficients so that the roots of the 
characteristic equation may be set at desired locations. 

In 1964, Siljak [2] extended Mitrovic's [1] , co- 
efficient plane method into the Parameter Plane method in 
which the two variable parameters may appear linearly in 
any or all coefficients of the characteristic equation. 

The parameter plane method provides the designer with a 
procedure for factoring a characteristic equation and 
displaying the results in a parameter plane diagram, and 
allows him to obtain information about system stability 
and how it is affected with parameter variations. 

The parameter plane method was further extended in 
1965 [4] to a general case in which the two variable 
parameters may appear nonlinearly in the coefficients of 
the characteristic equation. Another extension of the 
Parameter Plane method to include more than two variable 
parameters has been defined as the Parameter Space method 
[5] . 

When in general only two parameters (or coefficients) 
are involved, for any complex value of s or a pair of (£> w n ) 
values two simultaneous equations in two unknowns are ob- 
tained by setting the real and imaginary part of the 
characteristic equation equal to zero. 
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In 1967, Bowie [6] while working with the character- 
istic equation of an inertially stabilized vehicle on the 
parameter plane, found that a complete set of roots for 
the characteristic equation could not be determined in 
the real root area by using the existing parameter plane 
theory. He further investigated this problem and found 
that specific pairs of C , values, make the above men- 
tioned system of equations singular, which thus corresponds 
to a single straight line in the parameter plane, called 
a "singular line". 

The singular line added a new dimension to the para- 
meter plane because it provides the designer a way to hold 
a root fixed while varying others, by varying the two 
parameter values. Another important characteristic of 
singular lines is also the fact that they can be drawn in 
a parameter plane diagram by hand, without needing a 
computer, as in the case of drawing constant? or w n curves. 

Although an analytic study of the two previous men- 
tioned simultaneous equations, which are obtained by 
setting the real and imaginary part of the characteristic 

1 f 

equation equal to zero, makes obvious the existence of 
singular lines under certain conditions, they were not 
found earlier due to the standard computer programming 
practices and limitations. 

Singular line theory is based on the special case of 
linearly dependent parameter plane solution equations. 

From the standpoint of the design engineer it was 
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desirable that this theory be investigated in order to 
establish a standard general set of rules which would 
provide a method to obtain systems with singular lines. 

In this thesis this task is accomplished and possible 
applications of the singular line theory in control system 
design are investigated. 



17 



II. 



SINGULAR LINE THEORY 



A. SINGULAR CONDITIONS 



Consider the characteristic equation: 



n 



F(s) = 



E 




C 2 - 1) 



K=0 



where a ^ {k=0, 1, 2,...,n} , is a function of two variable 
parameters a and 3 , i.e., aj^gk (a, 3), and s is a point on 
the s plane which can be expressed in rectangular coordi- 
nates as: 



If for a given pair of w ns va l ues (where |C S | < 1 and 
w ns > 0 } j equation (2-1) has an infinite number of a and 3 real valued 
pairs which satisfy it, then for this equation singular conditions 
exist. In this case the point of the s plane which is defined by the 
pair £ , w ns is a singular point, and the characteristic equation is 
singular. 

The locus of all points (a, 3) on the Parameter Plane which 
satisfy equation (2-1) in the case of singular conditions is called 
the singular line. Each such line corresponds to a specific singular 
point or in other words to a specific pair °f £ s , w ns values. 




( 2 - 2 ) 
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B. SINGULAR LINE ANALYSIS IN THE LINEAR CASE 



A simple case, which was analyzed and which appears in 
many practical problems is the linear case. In this case 
the two variable parameters a and 3 appear in the coefficients 
of the characteristic polynomial in the general form: 

a K =3 K a+ c K 3 +d K (2-3) 

where b K , c K and d, c are real constants. 

1 . Analysis by Using Parameter Plane Equations 



Consider the characteristic equation (2-1), where 
a K is defined by equation (2-3). By introducing equation 
(2-2), into equation (2-1), and equating real and imaginary 
parts to zero, it can be shown [6] that: 



B. a + Cj g = 
B 2 a + c 2 6 = 

where 

»■=£; c-u Kb < “ n K v, to 

K = 0 

c -£ “ n K v. C!) 

K = 0 

Dx-f* (-l) K d a) K U x (?) 

1 i < n k- 1 

K = 0 



- Dx 

- d 2 




C 2 = 



t 

K = 0 



d 2 = 



L 



K = 0 



(- 1) K b oj K d U (?) 
Kn K J 

(- 1) K C u) K U (?) 
v J k n k 

(-i) K d oj K ;u (?) 

J k n k y J 



(2-4) 



(2-5) 
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The functions U k (^) and U k _ (£ ) are Chebyshev 
functions of the second kind, for which U o (S)-0, U-^(S)=1, 
and in general : 

u k C5)- 2CU k _ 1 CO - U k _ 2 CS) (2-6) 

The system C2-4) can be written in the Matrix form 

as : 









"" “ 


B 1 


C 1 


a 


' D 1 


B 2 


C 2 


6 


- D 2 



A x = B (2-7) 

When the solution of the system (2-7) is considered 
there are the following cases: 
a. No Solution 

When Rank (A) = 1 and Rank (A, B) = 2, which 
implies that: 




( 2 - 8 ) 



then the system has no solution, 
b. Unique Solution 

When Rank (A) = 2 and Rank (A, B) = 2 , which 
implies that: 




(2-9) 



then the system has a unique solution. 
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c. An Infinite Number of Solutions 



When Rank (A) = 1 and Rank (A, B) 

implies that 




1, which 



C2-10) 



then the system has an infinite number of solutions. There- 
fore in this case, since there are an infinite number of 
pairs of a, 3 values which satisfy the system (2-4) and in 
consequence the characteristic equation (2-1) , it implies 
that for this equation singular conditions exist. 



2 . Analysis by Using Complex Variable Theory 

The characteristic equation (2-1) is considered 
again. This equation can be rearranged in the form: 

f x (s) a + f 2 ( s ) 3 + f 3 (s) = 0 (2-11) 

where f^(s), f 2 (s) and f 3 (s) are polynomials of s. Intro- 
ducing equation (2-2) into equation (2-11), these poly- 
nomials will be in general a complex number with some 
magnitude and an argument. Thus they can be written as: 

f x (s) = r 1 e j0 l 

f 2 (s) = r 2 e j0 2 

f 3 (s) = r 3 ej 0 3 

Introducing the above equations into equation 
(2-11) yields: 

e-* 0 la + r£ e-* 0 2 3 + r 3 e^ 0 3 = 0 (2-12) 
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1 0 

Using the Euler's formula e J = cos0 + jsin0 , 
equation [2-12) becomes; / 

r 1 {cos0 1 +jsine 1 }oi+r 2 {cos0 2 + jsin0 2 }3+r 3 {cos0 3 +jsine 3 } - 0 

or 

{r^cos0^a+r 2 cos 0 2 3 + r 3 Cos 0 3 }+j{r^sin 0 ^a+r 2 sin0 2 3+r 3 sin0 3 } = 0 



Equating real and imaginary parts to zero yields 
the system: 



r^sin0^a+r 2 sin0 2 3=-r 3 sin0 3 

(2-13) 

r-^cos0^a+r 2 cos0 2 3=-r 3 cos0 3 

By applying the same reasoning as in paragraph one, 
the system (2-13) will have: 



a. No Solution 

If and only if: 



or 



r-^sin0^ 

r-jcose-^ 



r 2 sin0 2 , 
r 2 cos0 2 ' 



r^sin©^ 

r 3 COS0 3 



tan0^ = tan0 2 f tan0 3 

b. Unique Solution 

If and only if: 
tan0^ f tan0 2 j- tan0^ 

c. An Infinite Number of Solutions 
If and only if: 

tan©! = tan0 2 = tan0 3 



(2-14) 
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Equation (2-14) implies that for the character- 
istic equation (2-1) which has been rearranged in the 
form (2-11) singular conditions exist, if and only if: 



&1 = e 3 + k x n 

0 2 = 0 3 + k 2 n 
or 

/qcs) - /f,Cs) + kjH 
Acs ) - Ac*) * k 2 n 

where k^ and k 2 are integers. 



(2-15) 



3 . Singular Conditions for the Linear Case 

Summarizing previously derived results, the con- 
clusion is that, in the linear case where the character- 
istic equation is written in the general form: 



n 

F(s) = ]T a k s k = 0 

k=0 



where a^ = 3 ^ a + c k^ + ^k’ singular conditions exist if and 
only if: 




When the characteristic equation is considered 
in the general form: 

f 1 (s ) a + f 2 (s) B + f 3 (s) = 0 
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then singular conditions exist, if and only if; 
/f^s) -/4 3 Cs) ♦ k x n 

/f 2 Cs) -/ 3 (s) * k 2 n 



4 . Equation of a Singular Line 

When for a given pair of? s , w ns values singular 
conditions exist for the characteristic equation: 

n 



F( S ) -£ a k s k = 0 

k=0 



or 



f 1 (s)a + f 2 (s)3 + f 3 (s) = 0 

then, according to the previous analysis the two para- 
metersa and 3 will be related by one of the following 
equations, respectively: 



or 




(2-16) 



r T sin0_ r.sin0. r 7 cos0, r, cos0, 

3 = _ ± - _i ia = - -A £ - -1 i a ( 2 - 17 ) 

r 2 sin©2 r 2 sin9 2 r 2 cos© 2 r 2 cos0 2 



Each of these equations, is the equation of the 
singular line corresponding to the singular pair (? s , w ns ) 
which was considered. This singular line is a straight 
line on the parameter plane which is tangent to the constant 
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C s curve at the singular point [ 6 ]. It was found in this 
thesis, that this is also tangent to the constant w ns curve 
at the same point. 

5 . Methods to Test if a Point s is Singular 

When the characteristic equation is given and it 
is desired to determine if a given point ( £, w ) , ( | 5 | <1 
and w n > 0 ) , is singular the following methods can be used. 

a. Algebraic Methods 

By using the £ , values of the point under 
consideration, the coefficients of equations (2-5) can be 
evaluated. Then the given point will be singular if and 
only if equations ( 2 - 10 ) are satisfied. 

A second algebraic method consists of the 
evaluation of the arguments appearing in equation (2-15) , 
after writing the characteristic equation in the form of 
equation (2-11). Then the given point is singular if and 
only if equations (2-15) are satisfied. 

b. Graphical Method 

By this method the characteristic equation 
must be written in the form of equation ( 2 - 11 ) and the roots 
of the polynomials f-^(s), f 2 (s) and f^(s) must be evaluated. 
After these roots have been located on the s-plane, with 
the use of a spirule the arguments of the previous men- 
tioned polynomials for the point under consideration are 
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determined. Then this point will be singular if and only 
if, equations f2-15) are satisfied. 

c. Use of a Computer 

There are in general many ways by which a 
computer can be programmed to answer this problem. One 
such program is the << Singular Point Program >> which was 
developed with this thesis and which appears in Appendix 
H. In order for this program to be used in the case under 
consideration, the corresponding value of 5 of the point 
which is considered is used as an input data. Then at the 
computer print output the corresponding value of w n is 
checked if it is such that it makes the system (2-4) have 
an infinite number of solutions. When this is the case, 
the point ( C, w ) which was considered is a singular point, 
otherwise it is not singular. 

6 . Methods to Test if a Characteristic Equation is 

Singular 

In the previous paragraph a characteristic equation 
was considered to be given and it was desired to determine 
if a point of the s-plane was a singular point for this 
equation. In this paragraph the problem is reversed, i.e., 
a characteristic equation is considered to be given and it 
is desired to determine if there is at least one point of 
the s-plane which is singular. For this problem the fol- 
lowing methods can be applied: 
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a. Analytic Method 



By this method the characteristic equation is 
written in the form of equation (2-11). Then, the argu- 
ments of f^(s), f 2 (s) and f 3 (s ) are evaluated as functions 

of S and w , i . e . , 
n 7 

/f 1 Cs) = g-^ C,“> n ) 

/f 2 Cs) = g 2 c e,u> n ) (2-18) 

A 3 = ^ 3 ^ ?,u) n-* 



Introducing equations (2-18) 
(2-15) yields: 

g l(S,^ n ) = g 3 (C,o) n ) + k^n 

S 2^’%) = S 3^ ,a) n' ) + k 2 H 



into the system 



(2-19) 



When the system (2-19) is solved for the un- 
knowns C and u) n , any acceptable solution, i.e., any solution 

where the values of C and oj are real and also I £ I <1 

n 

and w > 0 , will define a pair of C, w values, which cor- 
n 7 c 7 n 

respond to a singular point of the characteristic equation. 
If at least one such point is found then the character- 
istic equation is singular. 

b. Graphical Method 

By this method the characteristic equation 
which is considered must be rearranged again in the form 
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of equation (2-11). Then the roots of the polynomials 
f^(s), f £ C s ) and fj(s) are evaluated and they are located 
on the corresponding points of the s-plane. By the use of 
a spirule a search on the s-plane must be performed for 
points which satisfy equations (2-15) . This search can be 
performed either by moving on constant C lines or constant 
w n cycles. This method will in general be useful when 
the polynomials of s which are involved in equation (2-11) 
are of low order and it is desired to find out if the given 
characteristic equation is singular for a restricted area 
of the s-plane which is of interest. 

c. Use of a Computer 

The Singular Point Program which was mentioned 
before and which was developed for this purpose corresponds 
to one among the many ways by which a computer can be pro- 
grammed to answer this problem. 
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III. 



SINGULAR LINE THEORY AND CONTROL SYSTEMS 



A. DEFINITIONS 



For a control system with two parameters a and g 
appearing in its characteristic equation the following 
definitions come as a consequence of the singular line 
theory which was previously developed. 

1 . Singular System 

A control system is singular, if and only if, 
its characteristic equation is singular. 

2 . Singular Characteristic Equation 

The characteristic equation of a control system 
is singular, if and only if, the system of the two 
simultaneous equations in two unknowns, a and (3, which is 
obtained by setting s = - jw \/l- into the character- 

istic equation and requiring the summation of reals and the 
summation of imaginaries to go to zero is singular, for at 
least one pair of ? > w n values, where | C | <1 and w n >0 . 
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3 . 



Singular Root 



A root of the characteristic equation of a 
control system is singular if it makes the characteristic 

equation singular. 

, / 

4 . Singular Point 



Each singular root corresponds to a specific 
point of the s plane. Such a point is called a singular 
point. 

B. SINGULAR LINE THEORY IN CONTROL SYSTEM DESIGN 

In general the analysis and design of any control 
system is based on an accepted model for the system under 
consideration. In the case of a linear control system of 
any order the no zero, two complex conjugate poles system 
model is used as a standard one. The accuracy of this 
approximation is tightly related to the dominant mode 
concept . 

At this point the use of the singular line theory 
seems to be of utmost importance since, if a system is 
singular the complex pair which provides the singularity 
remains fixed and the locations of the other poles may be 
adjusted to some extent because an infinite number of a 
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and 3 pairs is available for their adjustment. In addition 
some other performance criteria may be satisfied by an 
appropriate parameter variation. 

All these cases are analyzed and further investigated 
for other probable applications of the singular line theory 
in control system design. 

i 

C. SINGULAR LINE THEORY AND CONTROL SYSTEMS 
1 . Basic Requirements 



From all that has been stated in section II, it 
seems that in order for the singular line theory to be 
applied on control systems design and compensation the 
following basic requirements must be satisfied, i.e., the 
system under consideration must be singular, a singular 
root of the system characteristic equation must meet 
problem specifications when it is considered the dominant 
root, this root must be made dominant by adjusting the 
parameter a and 3 values, and finally the two parameters 
a and 3 must be adjustable. 

. Therefore, the initial problem in general will be to 
select the structure of the system or the required 
compensator for a given system, in order for the above 
basic requirements to be met. 

In the following paragraph a preliminary analysis of 
the form of the characteristic equation of a control system 
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was performed, as the first step in developing the re- 
quired structure of the control system to obtain the 
desired singular conditions. 

2 . The Characteristic Equation of a Control System 
in the Linear Case 



In general the characteristic equation of a con- 
trol system in the linear case will be in the form: 

F(s) = £ a k ° (3-1) 

k =0 



where a^ = b^. ct+c^. 3 +d^ and b^, c^. and d^ are real 
constants. Equation (3-1) can also be rearranged in the 
following general form: 

f 1 (s)a + f 2 ( s ) 3 +f 3 (s) = 0 (3-2) 

where f^(s), f 2 C s D and f 3 C s D are polynomials of s of any 
order. Assuming that f-^(s) and f 2 C s )) are not identical to 
zero, and that £ 3 ( 5 ) = 0 or f 3 (s) f 0 , yields the following 
two cases: 

a. £ 3 ( 5 ) ^ 0 

By assuming that £ 3 ( 5 ) f 0, dividing equation 
( 3 - 2 ) by £ 3 ( 5 ) yields: 



£ x (s) 

£ 3 (s) 



a 



f 2 Cs) 

f 3 (s) 



3+1 = 0 



(3-3) 



or 
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F i ( S ) q + F 2 (s) g + 1 = 0 



(3-4) 



where F-^s) and F 2 (s) are either a ratio of two polynomials 
of s or a polynomial of s. 

When the characteristic equation of a control 
system is considered to be in the form of equation (3-4), 
then the necessary and sufficient conditions for this 
equation to be singular, according to equations (2-15), 
are : 




where k^ and k 2 are integers. 

It must be noticed that according to the process 
which was followed in the derivation of equations (2-15), 
they are valid if either f-^(s), f 2 ( s D and fj( s ) are con- 
sidered polynomials of s or ratios of two polynomials of s. 
Based on this remark and the fact that the argument of 
unity is zero or an even multiple of n , equations (3-5) 
were derived from (2-15). 

Considering now the characteristic equation of a 
control system in the form of equation (3-4) the following 
cases may appear: 

(1) Case I: F^(s) f F 2 (s) # In this case a 

point s of the s-plane is singular for the charactereristic 
equation (3-4), if and only if, it does not lie on the real 
axis and satisfies equations (3-5). 
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In order to determine if a given point in such a 
case is singular or in order to find the singular points of 
the characteristic equation (3-4), if any, one of the 
methods which were described in the previous section can 
be used. In addition to these methods the following method 
which is related to the Root locus concept can be applied. 

The characteristic equation of the control system 
is : 

Fj.Cs) a + F 2 (s) 6 + 1 = 0 (3-6) 

Suppose that the Root loci of the following equations are 
drawn on the same drawing: 

A F 1 (s) = -1 

B F 2 (s) = -1 

C F 1 (s) = -1 

D F 2 (s) = -1 

where A and B are Root locus variable parameters which take 
values from zero to plus infinity and C and D are Root 
locus variable parameters which take values from zero to 
minus infinity. 

Based on the angle relation on which each Root 
locus is drawn, all the points of each of the above 
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Root loci will satisfy the following equations respec- 
tively: 

/F 1 (s ) = C2n ± 1) n 
/F 2 Cs ) - C2n ± 1) H 
/F 1 C,s ) = 2nn 
/f 2 (s ) = 2nn 

where n is an integer. Therefore, the singular points of 
the characteristic equation (3-6), i.e., the points of the 
s-plane which satisfy equations (3-5) simultaneously, will 
be the points where either the first or the third Root 
locus intersect with the second or the fourth Root locus. 

From the above described method it can be con- 
cluded that in the case under consideration either there 
are no singular points or there is a finite number of such 
points. Therefore, when in general F^(s) f F 2 (s) the 
characteristic equation of the control system will not 
always be singular and if it is singular the probability of 
finding a singular point which meets problem specifications 
as the operating point of the system will be very low. 

(2) Case II: F^(s) = F 2 (s). In this case 

the characteristic equation (3-4) of the control system 
becomes : 

F x (s) (o + 8) +1 = 0 (3-7) 
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08 ) 



Thus, the singular conditions given by (3-5) become: 



when there is only one angle relation to be satisfied, 
like that which is given by equation (3-8) , then there 
will be in general an infinite number of singular points. 
Considering equations (2-5) and (2-16), the equation of 
the singular line of any singular point will be in the 
general form: 



where C is a real constant. This implies that although 
in a singular case there will be an infinite number of 
pairs of a and 6 values which satisfy equation (3-7), the 
roots of this equation will be fixed for any of these 
pairs. For this reason this case can be considered as a 
trivial one under singular considerations and not useful 
for applications. 



/F 1 (s) = k x n 



As it will be shown in the following case, 



a + 6 = C 



(3) Case 111: F 2 (s) = F-^s). When F 2 (s) 



F-^ n (s), where n is an integer, then the characteristic 
equation (3-4) of the control system becomes: 



F 1 (s) a + F 1 n (s) 3 + 1 = 0 



(3-9) 



and the singular conditions given by (3-5) become: 




(3-10) 
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i.e., in this case a point s of the s-plane is singular, 
if and only if it does not lie on the real axis and satis- 
fies equation (3-10) . 

In order to determine if a point of the 
s-plane is singular for the characteristic equation (3-9) 
or in order to find the singular points of this equation 
one of the methods which were described in the previous 
section can be used. In addition to these, the Root locus 
method which was introduced in the first case can be used 
as follows: 

Suppose that the Root loci of the following 
equations are drawn: 

A F 1 (s) = -1 

B F 1 (s) = -1 

where A and B are Root locus variable parameters which 
take values from zero to plus infinity and from zero to 
minus infinity, respectively. 

Based on the angle relation on which each 
Root locus is drawn, all the points of each Root locus 
will satisfy the following angle relations respectively: 

/f 1 ( s) = (2k + l) n 
/F 1 ( s) = 2kn 

where k is an integer. Therefore, all the points of these 
two Root loci will also satisfy equation (3-10) . After 
these remarks it can be concluded that in the case under 
consideration there will be, in general, an infinite 
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number of singular points. Also the probability of finding 
a singular point which meets problem specifications as the 
operating point of the system will be relatively high. 

In the simplest case where n = 2, i.e., 

2 

F 2 CS) = (s) , the characteristic equation (3-9) becomes: 

F 1 (s) a+F 1 2 (s) 3+1 = 0 (3-11) 

and the necessary and sufficient conditions for this equa- 
tion to be singular will be given again by equation: 

ACs ) = kH (3-12) 

where k is an integer. 

(4) Case IV: F-^(s) = F 2 (s). in this case 

interchanging the two parameters a and 3 in the character- 
istic equation (3-4) yields the previous case. 

All the rest of the cases, which may appear can 
be considered to belong in Case I, and can be analyzed in 
a similar way. 

b. f 3 (s) = 0 

By assuming that f 3 (s) = 0, equation (3-2) 

becomes : 

f 1 (s) a+f 2 (s) 3-0 (3-13) 
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Since f^(s) and f 2 (s) are not equal to zero, dividing by 
£ 2 ( 5)8 yields: 



f x (s ) 

f^cl) 




1 = 0 



(3-14) 



or 

F, (s) - + 1 = 0 
1 3 



(3-15) 



The singular conditions for this case as it can 
be derived from equations (3-5) are given by the equation: 

/F^s) = kn (3-16) 

where k is a real integer or zero. It can be shown by the 
same procedure as was done in Case III in paragraph 2a(3), , 
that in this case there is in general an infinite number of 
singular points. Considering equation (2-5) and (2-16) the 
equation of the singular line of any singular point will be 
in the form: 



-f- = C (3-17) 

where C is a real constant. This implies that although in 
a singular case there will be an infinite number of pairs of 
a and 8 values which satisfy equation (3-15), the roots of 
this equation will be fixed for any of these pairs. For 
this reason this case also can be considered as a trivial 
one under singular considerations and not useful for 
applications . 
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3. 



Summary and Conclusions 



In the preliminary analysis of the character- 
istic equation of a control system, which was performed in 
the previous paragraph, it was found that there is a form 
of the characteristic equation, which in general has an 
infinite number of singular points. This is the following: 

F 1 (s) a+F 1 n (s) 6+1 = 0 (3-18) 

or in the simplest case where n = 2, the above character- 
istic equation becomes: 

F 1 (s) a+F 1 2 (s) 6+1 = 0 (3-19) 

where F^(s) f 0, and it is either a polynomial of s or a 
ratio of two polynomials of s. 

Therefore equation (3-18) corresponds to the most 
favorable case for the application of singular line theory 
in control system design and compensation. 

The next step is to find the structure of the 
system (or the required compensator for a given system) , 
which has a characteristic equation corresponding to the 
form of equation (3-18) or in the simplest case to equation 
(3-19) . This is done in the following section. 
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IV. 



ONE INPUT ONE OUTPUT LINEAR FEEDBACK SYSTEMS 



A. PRELIMINARY ANALYSIS 

In figure (4-1) the block diagram of a unity feed- 
back control system is shown. Suppose that it is desired 
to compensate the system by making use of singular line 
theory. Then according to the results of the previous 
section, the problem is to find the compensation structure 
or structures which provide a system with a characteristic 
equation in the general form: 

1 + a F (s) +3 F 2 (s) = 0 (4-1) 

i.e., a singular system. 

Considering a compensation structure then the adjustable 
parameters are in general all the compensator parameters, 
i.e., the poles, zeros and gain values, and the gain of the 
plant. Since the poles and zeros of a compensator can be 
either real or complex, a preliminary analysis is needed 
to choose the two parameters a and 6 from among the gain 
values . 

In the following sections which include a detailed 
analysis of the singular systems it is shown how one of 
these parameters, i.e., the parameter a, was introduced in 
the transfer function of a singular cascade compensator in 

such a way that its value affects its poles. ■ (See Section 
VI C and equation C6-18)). 
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In the preliminary analysis of the problem of finding 
the compensation structure or structures which provide a 
system with characteristic equation in the form of equation 
(4-1), i.e., a singular system, the five most commonly used 
compensation structures were analyzed as follows: 

1 . First Compensation Structure 

The first compensation structure which was con- 
sidered is that of figure (4-2) . 

The characteristic equation of this system is: 

1 + G (s) + G (s) H(s) = 0 (4-2) 

Comparison of equations (4-1) and (4-2) reveals 
that in order to obtain a singular system the following 
relations must be satisfied. 

G(s) = a Fg(s) 

H(s) = F gCs) (4 ' 3) 

i.e., the feedback compensator must have poles and zeros 
identical with those of the plant, and the two parameters 
a and 3 must be related with the gain values as is shown in 
equations (4-3). Introducing equations (4-3) into equation 
(4-2) yields: 

1 + a Fg (s) + 3 Fg 2 (s) = 0 (4-4) 

which is the characteristic equation of a singular system. 
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2 . 



Second Compensation Structure 



The second compensation structure which was 
considered is that of figure (4-3) . 

The characteristic equation of this system is: 

1 + G (s) G c Cs) = 0 (4-5) 

Comparison of equations (4-1) and (4-5) reveals 

that in this case there is no way by which a singular system 

can be obtained by defining the two parameters a and $ 

between the plant and cascade compensator gain values. 

In this case it was also noticed that when G(s) 

and G (s) are defined as: 
c 



and 



G (s) = kF a (s) 

o 



(s+Zj) (s+z 2 ) ... (s+z m ) 

s N (s+p 1 ) (s+p 2 ) ... (s+p n ) 



N 



G c (s) = ot+B kFg(s) 



a s (s+p 1 )...(s+p n )+0k(s+z 1 )...(s+z m ) 
s N (s+p : ) (s+p 2 ) . . . (s+pj 



the characteristic equation of the. system becomes: 

1 + akF (s)+Bk 2 F 2 (s) = 0 
S o 

Since the argument of the second and third term in 
the above equation are not affected by the appearance of the 
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2 

k and k factors in these terms respectively, the above 
equation is also singular. 

Based on the fact that the cascade compensator 
which was thus introduced in order for a singular system to 
be obtained has the same number of poles at the origin as 
the plant, this singular case (like the singular case, 
which was analyzed in Section V) has the basic disadvantage 
that the type number of the plant is not preserved, except 
in the case where the plant is type zero. 

3 . Third Compensation Structure 

The third compensation structure which was con- 
sidered is shown in figure (4-4) . 

Suppose that in general: 

G (s) = kFg (s) 

G C U) - k.F c (s) 

H(s) = k H F h (s) 

The characteristic equation of this system will 
be : 

1+G(s) G c (s) + G (s) G c (s) H(s) = 0 (4-6) 

Comparison of equations (4-6) and (4-1) reveals 
that in order to obtain a singular system the following 
relations must be satisfied: 

G ( s) G c (s) = kk c F g (s) F c (s) = aF g (s) FJs) (4-7) 

and 
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(4-8) 



H(s) = k H F h (s) = | F g Cs) F c Cs) 

Introducing equations (.4-7) and C4-8) into equa- 
tion (4-6) yields: 

l+cxF g (s) F c (s) + 8 F g 2 ( s ) F c 2 (s) = 0 (4-9) 

which is the characteristic equation of a singular system. 

This case can be considered (under singular consid- 
erations) identical with the first one of figure (4-2), 
when G(s) G c (s) is assumed to be one transfer function. 

4. Fourth Compensation Structure 

The fourth compensation structure which was con- 
sidered is that of figure (4-5). 

The characteristic equation of this system will be 
in general: 

1+G (s) H (s) + G (s) G c (s) = 0 (4-10) 

Comparison of equations (4-1) and (4-10) reveals 
that for a given plant i.e., given G(s) = kF CT (s), there is 

o 

an infinite number of ways by which G c (s) and H(s) can be 
defined in order for a singular system to be obtained. 

Based on the fact that certain types of feedback compensa- 
tion are in common use, e.g., velocity and acceleration 
feedback, only the following cases were selected to be 
analyzed as those which have the highest probability of 
application in practice. 



45 



a. HO) = a 



When HO) is a gain adjustment which is 
defined to be the parameter a, then in order for a sin- 
gular system to be obtained the following should be the 
transfer functions of figure (4-5) . 

GO) = kFg(s) 

G c (s) = eF g (s) (4-11) 

H (s) = a 

Introducing equations (4-11) into (4-10) 

yields : 

1+akFg (s) + ekF g 2 (s) = 0 (4-12) 

Because the factor k does not affect any argu- 
ment relation, equation (4-12) will correspond to a sin- 
gular system. 

b, H(s) - as 

When H(s) is defined to be a tachometer feed- 
back then in order for a singular system to be obtained 
the following should be the transfer functions of the 
blocks in figure (4-5). 

G(s) - kF g O) 

G c (s) - 6s 2 F g (s) (4-13) 

H(s) = as 
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Introducing equations (4-13) into equation 
(4-10) yields; 

1 + aks F ct (s) + 3ks 2 F 2 Cs) = 0 (4-14) 

which is the characteristic equation of a singular system. 

It should be noted that a basic steady state 
performance requirement of the system shown in figure (4-5) 
is that the cascade compensator G c (s) must transmit signal 
at zero frequency, i.e., its transfer function must not 
have any zero at the origin. But in this case it can hap- 
pen, if and only if, Fg(s) has at least two poles at the 
origin, which implies that the plant must be at least type 
two . 

c. HCs) = as 2 

When H(s) is defined to be an acceleration feed- 
back, then in order for a singular system to be obtained 
the following should be the transfer functions of the blocks 
shown in figure (4-5) . 

G (s) = kFg(s) 

G c (s) = as 4 FgCs) (4-15) 

H(s) = as 2 

Introducing equation (4-15) into equation (4-10) 

yields : 

1+aks 2 F g (s) + eks 4 F g 2 (s) = 0 C4-16) 

which is the characteristic equation of a singular system. 
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It should be also noted in this case that in 
order for the cascade compensator G c Cs) to transmit signal 
at zero frequency its transfer function must not have any 
zero at the origin, i.e., F (s) must have at least four 

O 

poles at the origin, which implies that the plant must be 
at least type four. Because such a plant usually cannot be 
found in practice this case will not be further considered. 

5 . Fifth Compensation Structure 



The last compensation structure which was con- 
sidered is that of figure (4-6) . 

By considering the product G(s) G^Cs) as one 
transfer function then under singular considerations this 
case turns to be identical with the previous one. 

B . SUMMARY 

In this section a preliminary analysis was performed 
on the most commonly used compensation structures in order 
to find which of these can provide a singular system, and 
under what conditions, i.e., a system with characteristic 
equation in the form of equation (4-1). 

The results of this analysis have indicated that from 
a theoretical point of view there are many ways by which a 
singular system can be obtained by properly selecting the 
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compensation structure and by properly defining the com- 
pensators and the two parameters a and 0 : 

The conclusion also from this analysis was that three 
singular cases must be considered for a further analysis, 
because they may be useful and applicable in control 
system design and compensation. These are the following: 
The first singular case corresponds to the compensation 
structure of figure (4-2) , where the feedback compensator 
and the two parameters a and 0 are defined by equation 
C4-3) . The second singular case corresponds to the com- 
pensation structure of figure (4-5) , where the compensators 
and the two parameters a and 0 are defined by equation 
(4-11) . The third singular case corresponds again to the 
compensation structure of figure (4-5), where the com- 
pensators and the two parameters a and 0 are defined by 
equation (4-13) , where the plant is considered to be at 
least type two. 

In the following three sections which refer to these 
three singular cases respectively, a detailed analysis is 
performed. 



49 



V. 



FIRST SINGULAR CASE ANALYSIS 



In the previous section the << first singular case >> 
was defined as that which corresponds to the compensation 
structure of figure (4-2) , where the associated transfer 
functions and the two parameters a and g are defined as 
follows : 



G (s) = aF g (s) = a 



(s+zj (S + Z 2 ) ... (s + z m ) 
s N (s+p 1 ) (s+p 2 ) ... (s+p n ) 



(5-1) 



and 



“(5) =§FgCs) -g 



(s + Zj) (s + z 2 ) 
s N (s+p 2 ) (s+p 2 ) 



(s + z ) 
v nr 



Cs+p n ) 



(5-2} 



A. SYSTEM OUTPUT AND ERROR ANALYSIS 



The system’s transfer function is: 

C (s) a F (s) 

= = e 

R(s) 1 +aF g (s) +3F g 2 (s) 



(5-3) 



or 



T(s) 



m n 

N a Q ( s+z j) Q (s+ Pk } 

at ji ^ xj n. m m 

s [^(s*P k ) +as [>♦* p^(s+Zj )+ep^(s+z_. f 
k=i k=i j =i j =i 



(5-4) 
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The output of the system is: 
C(s) = R[s) T (s) 



(5-5) 



and the error of the system is: 

E (s) = R(s) (1 - T (s) } (5-6) 

When the final value theorem was applied to equations 
(5-5) and (5-6) the steady state output and error were 
evaluated, assuming that all the poles of the transfer 
function T(s) were on the left half of the s-plane. 

The input signals considered are those, which are 
the most commonly used, i.e., step, ramp and parabolic 
input . 



1 . Plant Type Zero 



In this case where the plant is type zero, i.e., 
N 

N = 0, the factor s which appears in equation (5-4) be- 
comes unity. 



a. Step Input 

When a step input R(t) = Au(t) or R(s) = A/s 
was considered then: 



C(t) 

t+oc 



C (t) 
ss v 



lim s 
s->o 



m 



n 



-T(s) = 



An ^ 1 z . 1 1 

_ i=i Z j k= i 



n 



m 



1 1 p. 2 + a 1 1 z. 

k=l k j=i 3 






n m 

1 ip +3 "I Tz. 2 
k= i k j=i 3 
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and 



n 



m 



e (t) 

t^cc 



e (t) 
ss ^ ' 



= lim s- 

S + o 



(1-TCs) 



*n Pi! „ n 

}, V-l F k ■ l= 



i_Ll 



m n 



m 



~ip k 2+ ari2 n Pk+ 3n Zi 

k=i k j=i j k=i k j-i j 



b . Ramp Input 

When a ramp input RCt) = Btu(t) 
2 

or R(s) = B / s was considered then: 



C (t) = C ft) - lim s ~ T (s) = - 

t+oc S+O S Z 



and 



e(t) = e (t) = lim s -5, 



e ss Ct) = lim s — ^{l - T (s) } = « 
t - ** s^o s 



c. Parabolic Input 

2 

When a parabolic input R(t) = ct u(t) 

3 

or R(s) = 2c/s was considered then: 



C(t) = C (t) = lim s 4- T(s) = « 

t-^oc “ S+O S 



and 
e (t) 

t-^ 



e (t) 
ss v J 



= lim s 
s->-o 



2c 



(i-T(s) } = - 



2 . Plant Type One 

In this case where the plant was considered to 

be type one, i.e., N = 1, the steady error and output of 
the system were found as follows: 
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a. Step Input 



C(t) = C ft) = lim s £ T (s) = 0 

t+oc S+O * 

and 

e(t) = e (t) = lim s ^ {l-T(s)} = A 
t-*x 15 15 s+o 

b . Ramp Input 



C(t) = C (t) = lim s T(s) 

t^oc ^ S + O S Z 



m 



n 



n , n 



Ba Z . 



IZl i k= 1 k 

m 



n z. 



j = 1 



and 

e(t) = e ss (t) = lim s ~ {l-T(s)} = oc 

t-)“Oc S-^O 3 



c. Parabolic Input 



C(t) - C^(t) - 1 im s — j T(s) - “ 

s~*0 s 

and 

e(t) = e (t) = lim s {l-T(s)} = <* 
t-*-** s->-o s 



3 . Plant Type Two 



When a type two plant was considered, i.e., N = 2, 
the steady state error and output of the system were found 
as follows: 
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a. Step Input 



C (t) 
t-> « 

and 

e (t) 

t+cz 



C(t) 

t-> oc 

and 

e(t) 

t+“ 



C(t) 

t->oc 



and 

e(t) 

t->cc 



comes 

which 



C (t) = lira s - T(s) = 0 
ss s+o s 



e „(t) = lim s — {1-T (s) } = A 
s+o 



b . Ramp Input 



C (t) = lim s - T (s) = 0 

S->0 S' 



e ss (t) = lim s - 2 (l-T(s)} = 

s 



c. Parabolic Input 



C (t) = lim s — T (s) 

ss v ; 3 

S->0 s° 



m — 

. 2ca Pi z i Q 



m 

3 ri z . 

1=1 j 



: e ss^ t ^ ) = lim s ^ U-T(s)} = 



S->0 s 



Plant Type Three or Over Three 

In all such cases G s s C t ) = 0, and the error be- 

finite or infinity depending on the input signal 
is considered. 
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5. 



Remarks 



From the above results which are summarized in 
Table I of Appendix E, it can be seen that in general the 
type of the plant changes due to the form of the feedback 
path which is introduced in order for a singular system to 
be obtained. Therefore, compensation of a given system by 
this method yields a system which in general, does not 
respond to a given input as might be expected. 

Although some of the above cases which were 
analyzed may have a specific application in control sys- 
tem design and compensation, only the case where both 
output and error have finite values at steady state seems 
to be of a more general interest. 

The example which follows is referred to that 

case . 

B. EXAMPLE (5-1) 

Assume that in figure (4-2) the transfer functions of 
the plant and the feedback compensator are as follows: 

G(s) = aF (s) = a (5-7) 

S (s+i) (s + 2 ) (s + 4) 



and 



H(s) = | Fg (s) = £ 
v J a s v J a 



• S+3 

(s + l) (s+ 2 ) (s + *0 



(5-8) 
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then the characteristic equation of the compensated sys- 
tem is : 

1+ aF g (s) + 6 Fg 2 ( s ) = 0 (5-9) 

or after substituting equations (5-7) and (5-8) into (5-9) 
yields : 

s 6 +145 5 + ( 77+a) s 4 +(212 + 10a)s 3 +(308 +35a+B) s 2 + 

+ (224 + 50a + 66) s + (64 + 24a+9B) = 0 (5-10) 

By using the method of paragraph 6c, section II, the 
singular points of the system were found for several values 
of C. They are listed in Appendix A. One of these sin- 
gular points, i.e., the singular point (? s = 0.60, co ns = 
3.1104) was further considered. Then from equation (2-6) 
the Chebyshev functions of the second kind were evaluated 
and from equation (2-5) the coefficients B-^ , and 
were found to be: 

B 1 = -5.3085 
C ± = 0.6746 
= 41.7741 

Then from equation (2-16) the equation of the correspon- 
ding singular line was found to be: 

B = - 61.9254 + 7.8693a (5-11) 
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With a computer program then, the roots of the 
characteristic equation (5-10) were found for several 
different parameter pair values, which satisfy equation 
(5-11). They are listed in Appendix B. By examining 
these roots it can be seen that the singular root and a 
real root are fixed and the other three roots are varying 
for different pair values. 

In figure (5-1) the parameter plane diagram of this 
singular system is shown for some arbitrarily selected 
values of c, and 

In figure (5-2) the parameter plane diagram of the 
same singular system is shown for one constant £ and one 
constant w n curve, which correspond to the singular values 
of the previously considered singular point, i.e., = 

0.60 and w ns = 3.1104. On the same diagram the singular 
line of this singular point is also shown. It was noticed 
that these two curves and the singular line are tangent to 
each other on the same point of the parameter plane. 

For a further analysis of this system and of any 
other system which is referred to the << first singular 
case >>, the same procedure which is described in the next 
section and which is referred to the "second singular case", 
can be followed. 
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VI. 



SECOND SINGULAR CASE ANALYSIS 



In Section IV, the "second singular case" was defined 
as that which corresponds to the compensation structure of 
figure (6-1) where the associated transfer functions and 
the two parameters a and 3 are defined as follows: 



G(s) 



kF g (s) = k 



(s+Zj) (s + Z 2 ) . . . (s+Z^) 
s N (s+p) (s+p 2 ) ... (s+p n ) 



G c (s) 



SF g (s) = 3 



(s+Z 1 )(s+Z 2 ) . 
s N (s+p 1 ) (s+p 2 ) 



(s + Zj 
m 



. (s+p ) 



H(s) = a 



(6-1) 



( 6 - 2 ) 



(6-3) 



A. SYSTEM OUTPUT AND ERROR ANALYSIS 



The transfer function of the system is 



Tfsl - C(s) . G(s) G c Cs) 
1 tS ' ) RJJJ 



(6-4) 



1 + G (s) H(s) + G(s)G c (s) 



Introducing equations (6-1) through (6-3) 'into equa- 
tion (6-4) yields: 



m 

3k M (s + Z-) 2 

3 = 1 li 

?N n 2M n m m 2 

5 n (s + p k ) 2 + aks w n cs + p k )n(s + z j ) + 3 kn(s + z j ) 

k=i k=i j = i j = 1 

(6-5) 
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The output of the system is : 



CCs) = R(s) T (s) 



( 6 - 6 ) 



and the error of the system is: 



E(s) = R(.s) - C(s) = R(s) { 1 - T C.s) } 



(6-7) 



When the final value theorem was applied to equations 
(6-6) and (6-7) the steady state output and error values 
were evaluated, assuming that all the poles of the trans- 
fer function (6-5) lie on the left half plane. This 
analysis was based again on the three usually considered 
input signals, i.e., step, ramp and parabolic input. 

1 . Plant Type Zero 

In this case where the plant was type zero, i.e., 

N = 0, the steady state error and output of the system 
were as follows: 

a. Step Input 

When a step input R(t) = Au(t) or R(s) = A/s 
was considered then: m 



m 



C(t) = C fiS (t) = limsf T (s) 




2 




n 



n 



m 



m 
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and 



n 



e(t) 

t-** 



e (t) 
ss ' 




b . Ramp Input 



When a ramp input R(t) = Btu(t) 

2 

or R(s) = B/s was considered then: 

C(t) = C ft) = lim s — 2 T (s) = - 

t-*z S + O S 

and 

e (t) = e ft) = lim s \ (l-T(s)} = - 
t-** s+o s 



c. Parabolic Input 



2 

When a parabolic input R(t) = ct u(t) 



or R(s) 


= 2c / s 3 was 


considered then: 


C(t) = 


C (t) 
ss v J 


= lim s 
s->-o 


T(s) = - 
s 


and 








e(t) = 

t^oc 


e (t) 
ss v 


= lim s 
s^o 


{1-T (s) } = - 
s 
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2 . 



Plant Type One 



In this case where the plant was considered to 
be type one, i.e., N = 1, the steady state error and out 
put of the system were found to be as follows: 

a. Step Input 

When a step input was considered then: 

C(t) = C (t) = lim s £ T (s) = A 

t-*x S->0 



and 



e(t) 

t-w 



e ss (t) 



lim s - { 1-T (s) } = 0 
s+o ° 



b . Ramp Input 



When a ramp input was considered then: 



C(t) = C ft) = lim s \ T (s) = 
t+cc s s-*o s 



and 



e(t) 

t-^ 



e (t) 
ss 



u Ba O p 

lim s * 2 {l-T(s)}- 
s->° s 3 z j 



J 



= i 



c. Parabolic Input 



When a parabolic input was considered then 
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lim s 
s-^o 



cc 



Cft) - C ft) 



t->cc 



ss 




s 



and 




3 . Plant Type Two 

In this case where the plant was considered to 
be type two, i.e., N = 2 the steady state error and output 
of the system were evaluated as follows: 

a. Step Input 

When a step input was considered then: 



t+cc 



C (t) = C ft) 



lim s j T(s) = A 



and 



e(t) = e ss (t) = lim s | U-T(s)} = 0 




b . Ramp Input 



When a ramp input was considered then: 



C(t) = C„(t) = lim s \ T(s) 




cc 



and 



e(t) = e (t) = lim s U-T(s)} = 0 
t-w ss s+o s 
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c. Parabolic Input 



When a parabolic input was considered then: 



C(t) = e_„(t) = lim s T (s) 




oc 



and 



n 




e(t) = e (t) = lim s (l-T(s)} = 

"t* rr c 




4 . Remarks 

As a first step toward a further study of the 
"second singular case", an error and output analysis was 
performed. From these results which were obtained and 
which are summarized in Table II of Appendix E, it was 
noticed that in contradiction with the "first singular 
case" the type of the plant is preserved in this case, 
where the compensation structure and the associated com- 
pensators of figure (6-1) are introduced in order for a 
singular system to be obtained. 

■ f 

B. EXAMPLE (6-1) 

In this example a type one plant was considered, with 
transfer function. 



G (s) = kF (s) = 100 



1 



(6-8) 



s (s + 1) (s + 5) 
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Suppose that this plant is compensated by the com- 
pensation structure shown in figure C6-1) , then according 
to equations (6-2) and (6-3) the transfer functions of the 
compensators are: 

Gfs) = 3 F (s) = 3 ' ■ (6-9) 

s s (s + 1) (s + 5) 

and 

H(s) = a (6-10) 

The characteristic equation of the compensated system 
is : 

1 + G (s) H (s) + G (s) G (s) = 0 (6-11) 

Introducing equations (6-8) , (6-9) and (6-10) into equation 
(6-11) yields: 

s 6 + 12s 5 + 46s 4 +(60+100a)s 3 + (25+600a)s 2 

+500 as+1008 = 0 (6-12) 

Using the computer method of paragraph 6c, Section II, 
the singular points of the system were found for several 
values of z , . They are listed in Appendix C. 

One of these singular points was further considered, 
i.e., the singular point (? m 0.40, w ns = 0.97096). Then, 
corresponding to the damping ratio value, Chebyshev 
functions of the second kind were evaluated from equation 
(2-6) and from equations (2-5) and (2-16) the equation of 
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the singular line which corresponds to the singular point 
which was considered, was evaluated to be: 

3 = - 4.92497a + 0.25924 C6-13) 

Considering the last two terms of equation (6-12) a 
necessary condition for stability according to Routh’s 
first criterion is that a>0 and 3>0. Since a and 3 
parameters are related by equation (6-13) , it implies that 
a necessary condition for stability is: 

0< a < 0.05265 (6-14) 

By arbitrarily selecting a = 0.05 from the range of 
values defined by (6-14), the corresponding value of 3 
was found from equation (6-13) to be, 3 = 0.013. Then the 
roots of the characteristic equation (6-12) were evaluated, 
for this pair of parameter values, by using a computer 
program with the ZPOLR subroutine. These roots are: 

s 1 = - 0.3884 + j0.8899 s 2 = - 0.3884 - j0.8899 

s 2 = - 5.2141 s 4 = - 5.0132 

s 5 = - 0.9270 s 6 = - 0.0701 

i.e., one of the roots is the singular root (? = 0.40,w ns = 

0.97096) or s = - 0.3884 ± j0.8899. 
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In figure (6-2) the parameter plane diagram of this 



singular system is shown for some arbitrarily selected 
values of ? and w n . 

C. HARDWARE CONSIDERATIONS 

In figure (6-3) the block diagram of a unity feed- 
back system (plant) is shown. The open loop transfer 
function of this system is: 



Compensation of this system according to the "second 
singular case" scheme, yields the compensated singular 
system of figure (6-4). This compensation structure will 
be in general complex since it introduces a cascade and 
a feedback compensator. Therefore, a simplified version 
may also be desired. This modification is shown in figure 
(6-5) , where the feedback compensator (gain adjustment) 
has been taken from the plant and has been introduced at 
the cascade compensator G c (s) . Also the Root locus gain 
3 from G c (s) has been set as a gain adjustment in front of 
the minor loop. 

The characteristic equation of the systems in figure 
(6-5) and (6-6) are respectively: 



G(s) = kF g (s) 




(s+.zj ) (s.+ z 2 ). ( s + z m ) 



(6-15) 



1 + ak F (s) + 3k F 2 (s) = 0 

o o 



(6-16) 
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and 



1 + a F g (s) + (3k Fg 2 (s) = 0 (6-17) 

Comparison of these equations reveals that both are 
singular and they have the same singular points, but they 
will be in general different, except when k = 1, in which 
case they become identical. By this modification the 
whole compensation structure has been introduced in cas- 
cade with the plant, as is shown in figure (6-5). 

An equivalent system of figure (6-5) is shown in 
figure (6-6) where the whole compensation structure of 
figure (6-5) is replaced by an equivalent cascade com- 
pensator . 

The following considerations can now be stated con- 
cerning the compensation structures of figures (6-4) and 
(6-5) or (6-6) . First both are singular. Second with 
the compensation structure of figure (6-5) or (6-6) , the 
compensation problem has been simplified and better 
organized since it will be required to build a filter from 
distinct components or a filter in an integrated form 
which will be introduced in cascade with the plant. This 
is a simpler and in general a more realistic approach to 
the problem. In the case also where K = 1, i.e., the gain 
of the plant is equal to unity then these structures are 
also equivalent. When K f 1 then, although they will have 
the same singular points (because the singular conditions 



67 



to be satisfied are identical, in both cases), the 
corresponding equations of the singular lines will be 
different. This means that for the same singular point, 
different points of the parameter plane (a, 0) , will be 
associated and also the non-singular characteristic roots, 
will be in general different. 

The cascade compensator G gc (s) which was introduced 
in figure (6-6) has the transfer function: 



G 



sc 



(s) 



0 (S + Z j ) (s+ Z2-)' C s + Z m ) 



s (s+Pi)(s+p 2 ) ... (s+p n ) + a (s+z ! ) 



(6-18) 



(s + z J 
m 



where the subscripts s c stands for << Singular compen- 
sator >>. Assuming that a pair of a and 0 parameter pair 
values which satisfies problem specifications has been 
selected, then this compensator and its transfer function 
can be precisely determined. 

This modified singular compensation structure of 
figure (6-5) or (6-6) was further analyzed as the most 
applicable due to its simplicity under hardware considera- 
tions, in comparison with that of figure (6-4). It must 

* I 

be also noticed that in a specific application, where a 
gain feedback adjustment at the plant is possible or 
desired, e.g., when it is desired to eliminate the effect 
of some nonlinearities that the plant may have, then the 
initial compensation structure of figure (6-4) may be 
considered for application. 
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D. STABILITY ANALYSIS 



1 . Introduction 

One of the basic requirements of any control 
system is stability. The stability concept and related 
definitions can be found in reference [7]. The first 
question which arises when stability considerations are 
made for the singular structure of figure (6-6) is: What 

effect has the cascade singular compensator which has 
been introduced, on the stability of the initial uncom- 
pensated system? That is, in addition to the fact that 
this compensator makes the system singular, does it 
improve or impair the stability of the initial system? 

In either case the question which seeks an answer is: How 

is it possible from an initially stable or unstable 
system to obtain a stable singular system? The interpre- 
tations of curves and singular lines on a parameter plane 
diagram and determination of stability is not in general 
a simple process. 

Stability analysis by using the Root locus 
concept seems to be relatively easy and can give a more 
complete picture of the behavior of the system. The 
stability analysis as it can be performed by using the 
Root locus concept for the case of the compensated singular 
system is described below. 
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2 . 



Stability Analysis by the Root Locus Concept 



By substituting Fg (s) from equation (6-15) into 
equation (6-17) and after some manipulations the character- 
istic equation of the compensated singular system becomes: 



vr ^ 

a (s + z 1 )(s + z 2 ) . . . (s + z m ) s (s+p! ((s+p 2 ) . . . (s+p n ) +3k(s+z 1 ) . . 



2 N 



(s+ z m ) + s (s+p x ) (s+p 2 ) . ..(s+p n ) = 0 



(6-19) 



Assuming that the singular points of the singular 
compensated system have been found and that one of them 
has been selected as the operating point of the system, 
then the equation of the corresponding singular line can 
be evaluated from equations (2-5) and (2-16) . This 
equation will be in general in the form: 

3 = Aa+B 

where A and B are real constants. By substituting the 
parameter 3 from this equation (6-19) and after some 
manipulations, the characteristic equation of the singular 
compensated system becomes: 

2 2 2 

Ak(s+z x ) (s+z 2 ) ...(s+z m ) +s N (s+p 1 ) (s+p 2 ) . . . (s+P n ) (s+z 1 )(s+z 2 ) . . . (s+z m ) 

a. 2 2 2 2 N ~ 2 ~ (2 " “2 

Bk(s+z 1 ) (s+z 2 ) ...(s+z m ) +s (s+p x ) (s+p 2 ) ...(s+p n ) 

( 6 - 20 ) 
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Since A, B and k are known constants the numerator 
and denominator of the above equation can be factored, 
which in general yields: 



Cs + Zj' ) ( s + z 2 ) .... (s.+ z* ) 

a k- — = - 1 

(s+p{ ) (s+p 2 ' ) ... (s+p^ ) 



(6-21) 



where zl and p^f may be real or complex. By treating the 
parameter a as the Root locus variable parameter the Root 
locus of the compensated system can be drawn, either by 
the classical graphical methods or by the use of a com- 
puter. In the last case where a computer is used, usually 
it is possible to proceed directly from equation (6-19) to 
draw the Root locus of the compensated singular system. 

Since it is known that as the Root locus variable 
parameter a of equation (6-21) varies from zero to infinity 
the roots of the characteristic equation move on the Root 
locus segments, starting from the poles and terminating 
to zeros or at infinity when zeros are not available, and 
that the singular root and maybe one or more other roots 
of the characteristic equation (6-21) are fixed, indepen- 
dent of the variations of the parameter a, it implies that 
some zeros and poles of equation (6-21) have to be 
identical. Therefore, the numerator and denominator of 
equation (6-21) must consist of a number of common factors 
equal to the number of the fixed roots that the charac- 
teristic equation of the singular compensated system has. 
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This number is at least two, corresponding to the 
system's complex singular root. This fact makes the pro- 
cess of drawing the Root locus of the compensated 
singular system by the classical graphical methods, much 
simpler than it initially looks. Thus the Root locus of 
the singular compensated system consists of two or more 
fixed points of the s-plane plus the classical Root 
locus segments. 

Considering now system stability all the related 
questions can be answered from the Root locus study. For 
example, if any Root locus segment or isolated point lies 
entirely on the right half of the s-plane, then there is 
no pair of a and 3 parameter values for which the compen- 
sated system can be stable. If on the other hand all the 
Root locus segments lie on the left half of the s-plane 
as well as all the isolated points then all pairs of a 
and 3 parameter values will yield a stable system. Finally, 
if some Root locus segments cross the imaginary axis in one 
or more points , then the range of the parameter a values 
for which the compensated singular system is stable can be 
determined. The corresponding values of the parameter 3 
are then found from the equation of the singular line. 

The following example is an application of the above 
procedure and concepts which were introduced for the case 
of the singular compensated system. 
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3 . Example (6-2) 



In this example the uncompensated system (Plant) 
of figure (6r3), was considered to have transfer function: 



G (s) = kF (s) = k 

® s(s+l) (s+5) 



( 6 - 22 ) 



The characteristic equation of this system is: 



s 3 + 6s^ + 5s+k=0 



(6-23) 



The Root locus of this system is shown in figure (6-7). 

The gain at the stability limit is k=30. When this system 
is compensated by a singular cascade compensator the block 
diagram of the compensated system will be that which is 
shown in figure (6-5) or the equivalent of figure (6-6) . 
The singular cascade compensator which was thus introduced 
has a transfer function: 



G (s) = 
sc v J 



s 3 +6s^ + 5s+a 



(6-24) 



which was obtained from equation (6-18). 

The characteristic equation of the compensated 
singular system is: 



1 + G (s) G sc (s) = 0 



(6-25) 
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Introducing equations (6-22) and C 6 -24) into equation 
(6-25) yields: 

s^+12s^+46s^+(60+a)s^+(25+6a)s^+5as+k$ = 0 (6-26) 

Two values of the gain k were arbitrarily selected, 
i.e., k =20 and k =35, by which the uncompensated system 
becomes stable and unstable respectively. These cases 
were separately analyzed as follows: 

a. Uncompensated System Stable, k = 20 

For this value of k , the uncompensated 
system is stable and the characteristic equation (6-26) 
of the compensated singular system becomes: 



s 6 +12s 5 +46s 4 +(60+a)s 3 +(25+6a)s 2 +5as+20B = 0 



(6-27) 



In figure ( 6 - 8 ), the parameter plane diagram 
of the compensated singular system is shown for arbitrarily 
selected values of C and w n . By choosing s s =0.40 the 
corresponding value of w ns was found in the computer by using 
the << Singular Point Program>> as 1 ^ 5 = 0.97096. Therefore, 
a singular root of the compensated singular system is: 



s " " ?s “ns ± *> ns /l-s s = " 0.3884 ± j 0.8899 (6-28) 

From equations (2-5) and 2-16), the equation of the 
corresponding singular line was found to be: 

B = 0.24621a - 1.21242 (6-29) 

i.e., if the pair of a and 3 parameter values satisfies 

equation (6-29), then the characteristic equation (6-27) 
of the compensated singular system will have a fixed complex 
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root at the point s = - 0.3884 + j 0.8899. 

Considering now the coefficients of the zero 
and first power of s terms of the characteristic equation 
(6-27), a necessary condition for stability is that both 
these are greater than zero, which implies that <*> 0 and 
3>0. Based on this observation the roots of the character- 
istic equation (6-27) were found for a specific range of the 
parameter aand 3 values and with each parameter pair 
satisfying equation (6-29), by using the subroutine ZPOLR 
in a computer program. These roots are listed in Appendix 
D. From this list several general remarks can be made 
about the range of the a and 3 parameters in the stable 
region and the variation of the nonsingular roots in this 
region, i.e., the stable region in terms of the parameter 
a value is bounded between the values of 5 and 40 approx- 
imately, and that within this region the variations of 
the nonsingular roots is such that the dominance and 
sensitivity concepts are considerably affected, for dif- 
ferent parameter values. 

Substituting the value of 8 from equation (6-29) 
into equation (6-27) the characteristic equation of the 
compensated singular system becomes: 

s 6 +12s 5 +46s 4 +(60+a)s 3 +(25+6a)s 2 +5as+(4.9242a -24.2485) = 0 (6-30) 

Treating then the parameter a as the Root locus variable 
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parameter, equation (6-30) can be written as: 

s 3 + 6s 2 + 5s+ 4.9242 = - 1 (6-31) 

a — — 

s 6 +12s 5 +46s 4 +60s 3 +25s 2 -24.2485 
or by factoring numerator and denominator yields: 

(s+5. 2230) (s 2 +0. 7768s + 0.9428) 

a * * 1 1 

(s+5.2230) (s 2 +0 .7768s+0.9428) (s+4.7195) (s+1.8458) (s-0.5653) 

(6-32) 

The numerator and denominator of the left part 
of equation (6-32) have common factors corresponding to 
the real root s = -5.2230 and the complex singular root 
s = -0.3884 + j 0.8899. Therefore, the characteristic 
equation of the compensated singular system has a fixed 
real root and a fixed complex root at these points of 
the s-plane. 

The Root locus of this compensated singular 
system, which was obtained from equation (6-30). By the 
use of the R00TL0 subroutine in a computer program is 
shown in figure (6-9). 

At this point it must be noticed that, the 
remarks stated in paragraph two about the fixed roots 
that a singular system has, the analytical results shown 
in Appendix D, and the Root locus of figure (6-9) are 
in agreement. It must also be noticed that the fixed 
roots of the characteristic equation (6-30) are shown by 
dark cross marks in a circle, in figure (6-9). 
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b. Uncompensated System Unstable, k = 35 

When the gain of the uncompensated system 
(Plant) has the value of 35, then the uncompensated system 
is unstable since a pair of complex roots of the character- 
istic equation (6-23) lie on the right half of the s,-plane. 

In this case the characteristic equation (6-26) 
of the compensated singular system shown in figure (6-6) 
becomes : 

s^+12s^+46s^+ (60+a) s^+ (25+6a) s^+5as+35$ = 0 (6-33) 

In figure (6-10) the parameter plane diagram 
of this compensated singular system is shown. 

For the same value of C s as i n the previous 
case, i.e.,5 = 0.40, the corresponding value of w ns was 

found by using the "Singular Point Program" as w =0. 97096. 
Therefore, a singular root of the compensated singular 
system is s = -0.3884 - j 0.8899. The fact that the singular 
root remains the same as in the previous case, although the 
characteristic equation has been changed comes as a conse- 
quence of the singular line theory which was developed in 
section II and III, i.e., if the characteristic equation 
of the compensated singular system is considered to be in 
the form of equation (6-17), then in both cases Fg (s) is 
the same and only the value of k has been changed. Because 
the singular points according to the singular line theory 
are the points of the s-plane, which satisfy equation: 

/Fg (s ) = n II 
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where n is a real integer or zero, it implies that since 
Fg(s) is the same in both cases the singular points will 
remain the same. 

From equations (2-5) and (2-16) the equation of 
the corresponding singular line was found to be: 

3 = 0.14069a - 0.69281 (6-34) 

For the same values of the parameter ot, as those 

which were used in the previous case, the corresponding 

values of the parameter 8 were found from equation (6-34). 

For these formed pairs the roots of the characteristic 

equation (6-33) were evaluated by using the subroutine 

ZPOLR in a computer program. They are also listed in 

Appendix D. By comparing the roots of the characteristic 

equation (6-27) and (6-23), which are listed in Appendix D, 

it can be seen that they are identical for the same value 

of the parameter a. The reason for this fact is that when 

the singular point s = -0.3884 + j 0.8899 or (c = 0.40, 

s 

<o ns = 0.97096), is considered, the equation of the 
corresponding singular line as a function of k, according 
to equation (2-5) and (2-16) is: 

6 = i (4.9242a - 24.2485) (6-35) 

When the value of the parameter 8 from equation (6-35) is 
substituted in the general form of the characteristic 
equation of the compensated singular system, i.e., into 
equation (6-17) yields: 

1+a F g ( s ) + (4.9242a -24.2485) F g 2 (s) = 0 
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or after substituting F (s) from equation [ 6 - 22 ) and 

o 

making the necessary manipulations yields: 

s^+12s^+46s^+(60+a)s^+(25+6a)s^+5as+ 

(4.9242a - 24.2485) = 0 (6-36) 

i.e., the characteristic equation of the compensated 
singular system is independent of the value of k when the 
same singular point is considered. Therefore, for the 
same singular point and for the same value of the parameter 
a the compensated singular system will have the same 
characteristic roots independently of the value of k. As 
a consequence of these remarks not only for k = 20 but 
for any value of k >0 the Root locus of the compensated 
singular system will be identical with the Root locus 
shown in figure (6-9) when the singular point s = -0.3884 + 
j 0.8899 is considered. 

4. Summary and Remarks on the Stability of the 
Compensated Singular System 



When the plant to be compensated by a singular 
cascade compensator has no zeros in its transfer function 
then the parameter 8 and the Plant's gain k will appear 
only in the zero power of s term of the characteristic 
equation of the compensated singular system as can be 
seen from equation (6-19) . In addition, only the product 
of 8 and k will appear in this term, plus a constant when 
the plant is type zero. When a specific singular operating 
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point of the system is considered (S s > w ns )> then the para- 
meters a and 3 will be related by the equation of the 
singular line which corresponds to the above mentioned 
singular point. According to equation (2-5) and (2-16), 
the equation of the singular line will be in the form: 



where A and B are constants. From this equation it will 
be observed on the parameter plane that in the case under 
consideration when the value of k varies it will change 
the slope of the singular line which passes from the 
selected singular point. From this equation also it can be 
seen that the product 3k can be expressed as a function of 
the parameter a alone . By substituting this product in the 
characteristic equation (6-19) of the compensated singular 
system the whole equation will be expressed as a function 
of the parameter a. Thus the roots of the characteristic 
equation will depend on the value of a only. 



cascade compensator has one or more zeros in its transfer 
function then the parameter 3 and the plant's gain k will 
appear in three or more power of s terms of the character- 
istic equation of the compensated singular system as can 
be seen from equation (6-19). Because they will always 
appear together in a product form the coefficient of 
equation (2-5) can be written in the form: 



3 - | (A + Ba) 



(6-37) 



When the plant to be compensated by a singular 



n 




(6-38) 
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where C' is the constant coefficient of the 3k product, 

K 

in the j power of s term of equation (6-19). When a specific 
operating singular point of the compensated singular system 
is considered, i.e., C ,0J ns -* » t ^ ien " t ^ ie parameters a and 
3 will be related by the equation of the corresponding 
singular line of the above mentioned singular point. Then 
according to equation (2-5) and (2-16), the equation of the 
singular line will be: 

3 = ^ ( A * + B 'a) (.6-39) 

where A’ and B’ are constants. By continuing with the same 

reasoning as in the previous case, the same remarks were 

found to be applied here. 

Therefore it can be stated, that in general the 
stability of a compensated singular system when a specific 
singular point is considered, is independent of the value 
of the plant gain k and depends on the value of the para- 
meter a alone or in other words the stability of the compen- 

sated singular system can be interpreted in terms of the 
parameter a alone . The stable region also (if any) on the 
parameter plane will be bounded between Constanta lines. 

Thus when the compensated singular system is stable for 
some range of the parameter a values , it will also be 
stable for the same range of the parameter a values and 
for any value of the plant gain k. Therefore, when stability 
analysis is performed by the Root locus method where the 
parameter a is used as the Root locus variable parameter 

then the value of k can be set to any desired value as will 
be dictated from other design factors. 



81 



E. 



HOW TO OBTAIN A SINGULAR ROOT AT A DESIRED POINT 



OF THE S-PLANE 



1 . Introduction 

Compensation of a system by using a singular 
cascade compensator will be in general useful if the 
compensated system provides a singular root where the 
specifications of the problem dictate. Therefore, when 
designing by using the singular line theory the next 
problem which seeks an answer is: How a singular root can 
be obtained, where the system's specifications dictate? 

The general objective in all cases will be to 
make use of the advantage that, by varying the parameter 
values the singular root remains fixed and as other roots 
are varying the system's performance may be moved toward 
the system's specifications. 



2. Selection of the Singular Point 

Assume that a plant with transfer function: 



(s + z r ) (s + z 2 ) ... (s + z ) 

G(s) - kF (s) = k— n 2 

S (s+P! ) (s+P 2 ) ... (s+P n ) 



(6-40) 



is going to be compensated by a "singular cascade compensa- 
tor" whose transfer function according to equation (6-18) is 
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• • • 



(s + z ) (s + z ) ..... (s + Z ) 



q (s) = 3 i — — 

sc s N (s+P i ) (s+P 2 ) . . . (s+P n )+ a (s+z^) (s + z 2 ) ...(s + z m ) 

(6-41) 

then the characteristic equation of the compensated 
singular system will be: 



1 + aF (s) + Bk F 2 (s) = 0 (6-42) 

6 

According to equation (3-10) the singular points 
of the system will be the points of the s-plane which 
satisfy the following equation: 

/F g (s) = k x n (6-43) 

where k, is a real integer or zero. Therefore, the location 
of the singular points on the s-plane is determined only 
by the zeros and the poles of the plant, i.e., the singular 
points are a function of F (s) only, and in consequence 
they can be affected only by changing F CT (s). 

o 

In the case where among the singular points there 
is such a point which satisfies problem specifications as the 
operating point of the compensated singular system, then 
there is no need to change F (s) . Considering now the 

o 

case where such a singular point is not available then F (s) 
has to be changed in order to obtain a singular point where 
it is dictated by the problem's specifications. The fol- 
lowing examples were selected for illustration of the design 
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procedure in the above case, 
a. Example (6-3) 

Assume the plant to be compensated has the 
transfer function: 

G (s) = kF (s) = — — ^ (6-44) 

S s (s + 1) (s + 5) 

and that problem specifications, i.e., damping ratio, 
natural undamped frequency, settling time and peak 
overshoot require a dominant singular root at the 
point (c s = 0.4, w ns = 0.97096) or s s = - 0.3884 ± j0.8899. 

In figure (6-11) the poles of Fg (s) have been 
plotted on the s-plane as reference points, in order to 
measure the angles involved. The desired singular point 
is also located on the s-plane. By the use of a spirule 
it was found that: 

/F g (s) - - Uj + * 2 + <i> 3 } = - 180° 
therefore the point S s is actually a singular point since 
for this point equation (6-43) is satisfied for k^ = -1. 

This point also was found to be a singular 
point by the use of the "Singular point program" in 
example (6-1) . 

It must be noticed that another way to check 
if the point s ~ 0.4,^^=0.97096) is singular for the 
compensated singular system, is to draw the Root locus of 
the uncompensated system. Then if the desired operating 
point lies on the above Root locus (according to the 
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analysis which was developed in paragraph 2a(3), Section 
III), it is singular. Because this method is more 
laborious its use in general is not justified, except in the 
case where the desired operating point is not singular and 
some trade offs have to be done in order to decide if Fg(s) 
must be changed (by the appropriate initial compensator) 
or another operating point which lies on the Root locus of 
the uncompensated system can be selected as the operating 
point of the compensated singular system. 

b. Example (6-4) 

Assume now that the plant to be compensated 
has transfer function: 



and that problem specifications dictate for an operating 
singular point (C s =0 . 6 ,w ns =1) , or S s = -0.6±j 0.8. As in 
the previous example the poles and zeros of Fg(s) were plot- 
ted on the s plane as shown in figure (6-12). Then by the 
use of a spirule the argument of Fg(s) was evaluated and it 
was found that: 



Therefore the desired point is not singular since it does 
not satisfy equation (6-43) . This implies that Fg(s) must 
be changed somehow in order for the desired point to become 



G (s) = kF p.(s) = k — 

* s^ (s + 10) (s + 20) 



(6-45) 




(2xl27°+10.5°+5°) =- 247° 

(6-46) 
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singular. From a theoretical point of view it can be 
achieved as follows: 



(1) Introducing an Additional Zero. Suppose 
that an additional zero is introduced at the point s=-1.82. 
Then the angle contribution of this zero will be +67° as 
it can be seen in figure (6-13a) . Then, by defining the 
transfer function of the new obtained system after the 
above zero was introduced G z ( s ) , yields: 



which implies that the desired operating point (? s = 0.6, 
“ ns = 1) is singular for the system whose transfer function 



C2) Introducing an Additional Pole . Suppose 
that instead of a zero an additional pole is introduced 
at the point s=-0.58, then the angle contribution of this 
pole will be -113° as it is shown in figure (6-13b) . Then 
by defining the transfer function of the new system after 
the pole was introduced, Gp(s) yields: 



G z (s) = kF z (s) = k{ (s + 1 . 82) Fg (s) } 



and 





k (s + 5) (s+1 .82) 

s 2 (s+10) (s+20) 



G p (s) = k F p (s) = k { 




and 
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which implies that the desired operating point (£ s =0.6, 
to = 1) is singular for the system with transfer function: 

H.S , r* 

s + 5 



G (s) = k 

P s (s + 0. 58) (s + 10) (s + 20) 

( 3 ) Introducing an Additional Pole and Zero . 

Suppose now that an additional pole and a zero are 
introduced at the transfer function of the initially 
given plant. This pair of points can be chosen in such 
a way that their net angle contribution will make 
the desired singular point ( C = 0.60,w ns = l) actually 
singular. For example if the pole is chosen at the point 
s = - 16 and the zero at the point s = - 1.66 as it is shown 
in figure (6-13c), then by defining the transfer function of 
the new system after this pole and zero have been introduced, 
G zp (s) yields: 

G zpCs) - k F zd Cs) = k { s ^; 66 F g (s)> 



and 



zp 



(s) = / F (s) + 73 



As 7 



g 



s + 16 

° - 6 ° = 



180‘ 



which implies that the desired operating point is singular 

for the system with transfer function: 

(s + 5) (s+1.66) 

G (s) = k 

zp 

s 2 (s + 10) (s + 16) (s + 20 
3 . Summary and Remarks 

Considering in general an uncompensated plant with 
open loop transfer function G(s) = k Fg(s), where Fg(s) 
is a ratio of any two polynomials of s, then the argument 
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of F (s) evaluated at any point of the s-plane will always 

§ 

lie between the limits: 

k x n < / Fg ( s) < C.k-^+1) n 

where k^ can be any real integer or zero. Therefore, the 
maximum angle correction which Fg(s) may require in order 
to satisfy equation (6-43) for a desired point of the 
s-plane will be less than 180° , which implies that a 
single zero or a single pole or a zero and a pole or finally 
any number of zeros and poles can be properly selected and 
introduced to Fg ( s ) in order to make the argument of Fg(s) 
satisfying the angle requirement for singular conditions, 
for a specific singular point, which is desired to be the 
operating point of the compensated singular system. 

These poles and zeros in general may be real or 
complex. In consequence it implies that when such a change 
of Fg(s) is required an initial cascade compensator must be 
introduced, before a "singular cascade compensator" is 
introduced. 

Since a compensator with only zeros actually can 
not be built, the case of introducing one or more zeros 
alone has to be ignored as unrealistic. Therefore, among 
the rest of the possibilities available an optimal selection 
of the initial cascade compensator has to be done. This 
selection can be based on noise considerations, problem 
specifications and the "dominance of the singular root” 
concept which is described in the next part of this section, 
i.e., except the fact that the Root locus of the uncompensated 
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system will be forced to pass through the point where the 
desired singular root is located, the shape also of the 
resulting Root locus must be such that it promises some 
range of dominancy of the singular root. (See Example 
( 6 . 6 )). 

It must also be noticed that as it will be shown 
in the following analysis when an initial cascade compensa- 
tor is required, then this requirement will actually affect 
only the transfer function of the singular cascade compensa- 
tor and neither the complexity of the rest of the design 
problem nor the compensation structure (See Example (6-6) 
and figure ( 6 - 28 c)). Under these considerations an initial 
cascade compensator with only zero's may also be considered. 

F. DOMINANCE OF THE SINGULAR ROOT 
1 . Dominant Root Concept 

The concept of the dominant root has wide applica- 
tions in analysis and design of linear systems. Considering 
a system with a characteristic equation of order n, then 
there will be n characteristic roots for this system. If 
one of these roots dominates the transient response, then 
this root will be called dominant. When the dominant root 
of a system is a complex root, s = - j w n /l - 5 , then the 
overall response of the system for a given input may be 
closely approximated from the corresponding standard second 

order system graphs (i.e., the second order system whose 

2 2 

the characteristic equation is s + 2£co s+w =0 and has no 

n n n 

zeros) . 
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2 . Dominance of the Singular Root 



When a singular root has been determined from 
the problem specifications, then this singular root may 
be available from the given plant or not. In the last 
case an initial cascade compensator will be required, in 
order for the desired singular root to become available. 

The design process of this initial cascade compensator 
was previously illustrated. 

In either case it will also be desired for the 
above considered singular root to be dominant for one or 
more (depends from the specific application) parameter pair 
values. The design process in all possible cases toward 
the dominancy of the singular root will be illustrated by 
the following examples: 

3 . Example (6-5) 

In this example the Plant of example (6-2) was 
again considered with the Root locus gain k = 20, i.e.: 

G(s) = k F g (s) = 20 (6-47) 

Since the value of the plant's gain does not 
affect the location of the singular roots on the s-plane, 
as was shown in example (6-2), the plant under consideration 
when it is compensated by a "singular cascade compensator" 
will have the same singular roots with these which are 
listed in Appendix C for the example (6-1) and (6-2). Some 
of these roots were considered and their dominancy was 
investigated. 
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In figure (6-14a) and (6-14b) the compensated 
singular system is shown, where in figure C6-14a) the 
singular compensator is considered to be in integrated 
form, and in figure (6-14b) the same singular compensator 
is considered to consist of distinct components. 

The singular cascade compensator which thus has 
been introduced according to equation (6-18) has transfer 
function: 

G rsl - B ' ' C6-48) 

sc lSJ s(s + l) (s + 5) + a 

The characteristic equation of the compensated 
singular system was found to be: 

s^+12s^ + 46s^+ (60+ct) s^+ (25 + 6a) s^ + 5as +208 = 0 (6-49) 

a. Singular Root (C = 0-40, w ns = 0.97096) 

When the singular root (? s = B.40,u ns = 0.97096) 
was considered, then the problem which was investigated was: 
Are there any pairs of parameter a and 8 values for which this 
singular root becomes dominant? 

From equations (2-5) and (2-16) the equation 
of the corresponding singular line was found to be: 

8 * 0.24631a - 1.21242 6-50) 

Introducing equation (6-50) into equation 
(6-49) the characteristic equation of the compensated 
singular system in terms of the parameter a becomes: 



91 



(6-51) 



s^+12s^+46s^+(60+a)s^+(25+6a)s^+5as 
+ (4. 9242a - 24.2485) = 0 

Treating then the parameter a as the Root locus 
variable parameter, the Root locus of the compensated 
singular system for the singular point under consideration 
was drawn in the computer. It is shown in figure (6-9) . 

The compensated singular system has three fixed roots, i.e., 
the singular root s s =-0.3884 ±j 0.8899 and a real root 
s=-5.2230, which are shown with dark circle marks at the 
corresponding locations of figure (6-9) . The compensated 
singular system also has three other roots whose location 
varies on the corresponding Root locus segments according 
to the value of the parameter a . 

Examining the roots of the characteristic equation 
C6-51) in relation with the Root locus of figure (6-9) 
it can be seen that the value of the plant gain does not 
affect the location of the roots of the characteristic 
equation of the compensated singular system and that the 
location of the three varying roots is affected only by 
the value of the parameter a . The dominancy also of the 
singular root can be affected only by the other two 
characteristic roots which lie on the same Root locus seg- 
ments with the singular root. These roots are either real 
or complex according to the value of the parameter a . 
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It was found that in the stable region the value 
of the parameter a is approximately betwen the values of 
five and fourty. When a = 9.85 then the two varying roots 
in the vicinity of the singular root become identical with 
the singular root, i.e., for the value of the parameter 
a = 9.85 the singular root is repeated. It was also found 
that a=9.85 not only the singular root but and the real 
fixed root ats=-5.2230 was also repeated. 

In an attempt to find the values of the parameter 
a for which the singular root is dominant Cor in other 
words the system’s performance can be well predicted from 
the corresponding to the singular root second order model) 
several runs with different parameter pair values each, 
were performed in the computer using a DSL (Digital 
Simulation Language) program. In this program the compensa- 
ted singular system shown in figure (6-14b) was implemented. 

When the second order model corresponding to the 

singular root (C = 0.40,u) =0.97096 was considered with 

b s ’ ns 

transfer function: 



C(s) = 
R(s) 



0) 



ns 



0.9428 



(6-52) 



S +2C s u ns S + “ns 2 s 2 +0 . 77685S+0.9428 



its transient response characteristic values were calculated. 
These values together with those which were obtained from 
the computer for the compensated singular system are 
listed in Table I of Appendix F. 
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By examining this table the following remarks can be stated: 



Cl) There is no pair of parameter values for 
which all the transient response characteristics of the 
compensated singular system can be well predicted from 
the values corresponding to the singular root second order 
model . 

C2) For all the parameter pairs both the 
singular root and the other root in the vicinity of the 
singular root contribute to the transient response of the 
compensated singular system. 

(3) Considering that in practice the only 
transient response characteristic values for which the 
designer is in general interested, are the settling time 
and the percent overshoot, then from this point of view, 
it can be said that the singular root is dominant for 
a =7.8 and g = 0.708. Assuming also that there is some 
tolerance in the system performance then these two parameter 
pair values can vary within some range according to the 
tolerance values. 

The conclusion from this example was 
that when dominant mode design is used for compensation of 
the given plant then the singular line theory provides in 
general an acceptable solution. In figure (6-15) the time 
response of the compensated singular system is shown for 
a unit step input and for the two parameters having the 
values of a= 7.8 and 3 = 0.71. In figure (6-16) the time 
response of the compensated singular system is shown for 
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the parameter pair values a= 9.85 and 8 = 1.21 for which the 
singular root and the fixed real root are repeated. This 
fact as can be seen from this figure causes a very high 
overshoot value. In figure (6-17) the time response of 
the compensated singular system is shown for such a pair 
of values, i.e.,a =20 and 8= 3.71, for which the contribution 
of the singular root and the other root in the vicinity of 
the singular root can be recognized. In figure (6-18) the 
time response of the compensated singular system is shown 
for the case where the characteristic root in the vicinity 
of the singular root lies on the real axis. In this case 
because this root is closer to the origin and to the imaginary 
axis than the singular root, it is dominant and the transient 
response of the compensated singular system has an exponential 
form. 

b. Singular Root (?_ = 0.30,w =1.15268) 

5 1 1 o 

The second singular root from the table of 
Appendix C which was examined for dominancy was the singular 
root (5 S = 0.30,w ns = 1.15268) or s =-0.34580 ± j 1.09959. 

From equations (2-5) and (2-16) the equation 
of the corresponding singular line was found to be: 

8 = 0.35266a - 2.48732 (6-53) 

Introducing equation (6-53) into equation 
(6-49) the characteristic equation of the compensated 
singular system in terms of the parameter abecomes: 

s 6 +12s 5 +46s 4 (60+a)s 3 +(25+6a)s 2 +5as+7. 0532a - 49.7464 = 0 

(6-54) 
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Treating then the parametera as the Root locus 
variable parameter, the Root locus of the compensated 
singular system for the singular root under consideration 
was drawn in the computer. It is shown in figure (6-19) . 
This Root locus as can be seen from this figure consists 
of three points of the s-plane where the three fixed roots 
of the compensated singular system lie and the classical 
Root locus segments. The above mentioned three fixed roots 
are shown with dark circle marks at the corresponding 
locations , i . e . , 

s x = - 5.3084 
s 2 = - 0.34580+ jl. 09959 
s 3 = - 0.34580 -jl. 09959 

In an attempt to find the values of the parameter 

a for which the singular root is dominant several runs with 

different parameter pair values each, were performed in the 

computer using the DSL simulation program. The transient 

response characteristic values which were obtained from 

the computer, together with these corresponding to the 

singular root (t =0.30,w = 1.15268) second order model 

s ns 

are listed in Table II of Appendix F. 

By examining this table, the same general remarks 
and conclusions were derived in this case as in the previous 
one. The basic remark again was that there is a pair of 
parameter values, i.e., a= 10.6, 8= 1.25, for which the 
singular root is dominant from a practical point of view. 

In figure (6-20) the time response of the compensated 
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singular system is shown for this pair of parameter values. 
In figure also ( 6 - 21 ) through (6-23), the time response of 
the compensated singular system is shown for several 
characteristic pairs of parameter values. 

c .' Singular Root (C = 0.75,to ns = 0.60667) 

The third singular root which was examined 
for dominancy was (£ = 0.75 , 0J nS = 0.60667) or 

s = -0.45500 ± jO. 40127. 

From equations (2-5) and (2-16) the equation 
of the corresponding singular line was found to be: 

3 = .0.0936 70 - 0.17547 (6-55) 

Introducing equation (6-55) into equation 
(6-49) the characteristic equation of the compensated 
singular system in terms of the parameter a becomes : 

s 6 + 12s 5 + 46s 4 +(60 + a)s 3 +(25 + 6a)s 2 + 5as + 1.8734a-3.5094 = 0 (6-56) 

Treating then the parameter aas the Root locus 
variable parameter, the Root locus of the compensated 
singular system for the singular root under consideration 
was drawn in the computer. It is shown in figure (6-24). 
This Root locus (as can be seen from figure (6-24) also 
consists of three points of the s-plane where the three 
fixed roots of the compensated singular system lie and the 
classical Root locus segments. The three points where the 
three fixed roots are located are shown with dark circle 
marks at the corresponding locations, i.e.: 
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5.0899 



S 1 

s 2 = - 0.45500+ jO. 40127 
s 3 = - 0.45500- jO. 40127 

By the same procedure used in the previous cases 
the same general remarks and conclusions were derived in 
this case also. The basic remark again was that there is 
a pair of parameter values a and 8 for which the singular 
root is dominant under a practical point of view. This pair 
as can be seen from the Table III of Appendix F is: ct = 3.6, 

3 = 0.162. In this table the transient response character- 
istic values which were obtained from the computer together 
with these corresponding to the singular root (c g = 0.75, 
w = 0.60667) second order model are listed. In figure 
(6-25) the time response of the compensated singular sys- 
tem is shown for the previous mentioned parameter pair 
values . 

d. Remarks 

It must be noted here that the proceeding 
design was available for the particular system. It may 
not be generally available or may not provide an accept- 
able solution to some problems. 

4. Root Loci Correlation 



In the previous example if the desired singular 
root of the compensated singular system could not be made 
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dominant by a pair of parameter values, then an initial 
cascade compensator would be required, as has been already 
stated in paragraph two. 

Since the initial cascade compensator is applied 
to the uncompensated system, i.e., the system before the 
singular cascade compensator has been introduced, in order 
to affect the compensated singular system, i.e., the 
system after the singular cascade compensator has been 
introduced, a prior knowledge of this effect is necessary 
in order to avoid an undetermined number of trials. 

This prior required knowledge or the criterion 
for the selection of the initial cascade compensator can 
be based on the correlation between the Root locus of the 
uncompensated system and the Root locus of the compensated 
singular system, as can be seen from the following analysis. 

When an uncompensated linear control system with 
unity feedback is considered, its characteristic equation 
will in general have the form: 



(s + z^ (s + z -) . . . (s+z m ) 
s N (s+p 1 ) (s+p 2 ) ... (s+p n ) 



k F g Cs) = - 1 



( 6 - 57 ) 



where k Fg(s) is the open loop transfer function and k 
is the Root locus gain. The Root locus of this uncompen- 
sated system which is drawn with the Root locus gain k as 

the Root locus variable parameter consists of all the 
points of the s-plane which satisfy the angle requirement: 
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/\ Cs) = ( 2 n± 1 ) II 

Therefore any point of the Root locus of the uncom- 

pensated system as it can be seen from equation (6-57) 
and the previous angle relation satisfies the following 
two relations: 



and 



or 



and 



k F g (S RL ) = - 1 



/V s rl) - 



F g( S RL) < 



/Fg(S RL ) * (2 ni U n 



(6-58) 



(6-59) 



where n is a real integer. 

When the system which was considered is com- 
pensated by a singular cascade compensator the charac- 
teristic equation of the compensated singular system 
according to equation (6-17) will be: 



l+o F g (s) +$ k F g ( s ) = 0 



(6-60) 



Then according to equation (3-12) a point of 
s-plane is a singular point for the compensated singular 
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system if and only if it does not lie on the real axis 
and satisfies equation: 



/F g fs) = nk 



C6-61) 



where n is a real integer or zero. Therefore all the 
points of the Root locus of the uncompensated system, which 
do not lie on the real axis (this restriction is implied 
from the definition of the singular point in Section II, 

A) are singular points because they satisfy equation (6-59) 
and in consequence equation (6-61) . 



point (C s > w ns ) is selected as the operating point for the 
compensated singular system. For this point equation 
(6-58) becomes: 



(i.e., the equation which defines the relation between the 
two parameters a and 8 for which the selected singular 
point will always be a root of the characteristic equation) 
can be found either by using equations (2-5) and (2-16) 
or by direct substitution of the selected singular point 
in equation (6-60), which yields: 



Suppose that one of these points, i.e., the 




0 



(6-62) 



The equation of the corresponding singular line 



1+a k F 8 2 CC s .“ ns J = 0 



(6-63) 



101 



Solving this equation for 3k and substituting 



it in equation (6-60) yields: 



1+a F C.r ) 

l+« F g (.s) - s__ns F g Z (s) = 0 

S v s ’ ns' 



C6-64) 



Equation (6-64) is the characteristic equation 
of the compensated singular system as a function of the 
parameter a, when the singular point is considered 

to be the operating point of the compensated singular sys- 
tem. This equation can be written in the form: 



a 



tF g^s' M ns J - V s ) > 
< F g< 5 s*"ns 3 - ^ 



(6-65) 



( 6 - 66 ) 

F ct (C ,oj )(s + z )(s + z. ) . . . ( s + z ) 

g ^ s ns' v i M 2 m' 

F g (C ,0) )s N (s+p )(s+p ) • . . (s+p ) + (s + z ) (s + z )...(s+z ) 

o b lib 12 il 1 2 li 



Considering the Root locus of the uncompensated 
system, the Root locus of the compensated singular system 
and the corresponding equations, the following remarks 
and conclusions can be stated. 
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a. Number of Fixed Roots 



From equation (6-65) , the common factor in 
the numerator and the denominator on the left side of 
equation, implies that the compensated singular system 
has a number of fixed roots equal to the order of the 
characteristic equation of the uncompensated system, e.g., 
when the order of the characteristic equation of the 
uncompensated system is three, there will be three fixed 
roots, two of which have to be the complex pair of the 
singular roots. 

b. Root Loci Correlation 

The angle relation which the points of the 
Root locus of the compensated singular system satisfy is: 



/ Fg U s >a ns )F g (S) 

j F g ^s >a) ns^ +F g ^ 



( 2n±l) n 



(6-67) 



Then according to equations (6-58) and (6-62) all the 
points of the Root locus of the uncompensated system satis- 
fy equation (6-67) . Therefore all the points of the Root 
locus of the uncompensated system will be also points of 
the Root locus of the compensated singular system. 

Second, from equation (6-57) and (6-66) it 
is observed that they have, at the expression on the left 
side, identical zeros. They also have the same number of 
poles . 
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The conclusion from the above remarks is 
that the Root locus of the compensated singular system will 
be identical with the Root locus of the uncompensated sys- 
tem with the only difference that the Root locus segments 
of the compensated singular system may be further extended 
from the pole side end points, where the corresponding 
Root locus segments of the uncompensated system terminate, 
until they reach the poles of the compensated singular 
system . 

Examining the figures (6-7), (6-9), (6-21) 
and (6-25) they confirm all the above derived remarks and 
conclusions. Therefore, if an initial cascade compensator 
is introduced to the uncompensated system this will affect 
the Root locus of the uncompensated system and the Root 
locus of the compensated singular system in the same way, 
i.e., the prior knowledge of the effect of the initial 
compensator on the compensated singular system can be 
obtained by studying the effect of the initial compensator 
on the uncompensated system. 

5 . Range of Dominancy of the Singular Root 

When a singular compensated system is considered 
and an operating singular point, then the parameters a and 
3 will be related by the equation of the corresponding 
singular line. In general there will be two ranges of 
a and 3 parameter values of primary interest. The first 
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one will be that for which the compensated singular system 
is stable and the second one will be that within the pre- 
vious one, in which the singular root is dominant. In 
general the wider these ranges the more desirable will be 
the solution. Therefore when an initial cascade compensator 
will be introduced to the uncompensated system in order to 
make the singular root dominant, it must be selected in 
such a way that it will also give a wide range of dominancy 
of the singular root. 

It must be realized that in a self-adaptive 
singular system where the two parameters a and 3 are 
automatically adjusted so that the singular point will 
remain on the singular line during the operation of the 
system, the range of dominancy of the singular root will 
be very important. 

The following example is mainly referred to the 
range of dominancy of the singular root. 

6 . Example (6-6) 

In this example the same plant which is shown 
in figure (6-26) was again considered. Its transfer 
function is: 

GCs) = s"( s + 1 y ( s"+" 5") 
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For this example, it was assumed that prohlem specifica- 
tions require an operating point for the system, i.e., an 
operating dominant singular root, with damping ratio 
S = 0.60, and natural frequency = 1.2 rad/sec. A wide 
range of dominancy was also assumed to be required for the 
above singular root. 

The Root locus of the plant with the Root locus 
gain used as the Root locus variable parameter is shown in 
figure (6-29), where the location of the desired singular 
root is also shown. By studying this figure and taking 
into account the remarks which were stated in part VI E 
and VI F, the conclusion was that an initial cascade com- 
pensator is needed for the following • two reasons: The 

first is that since the desired singular root does not 
lie on the Root locus of the plant this Root locus must be 
reshaped in order to pass through the point where the 
desired singular root is located, and the second reason is 
that the shape of the Root locus must be changed in order 
for this singular root to become dominant with a wide range 
of dominancy, as it is required from the problem specifi- 
cations . 

After the above remarks and taking into account 
the noise problem it was decided to compensate the plant 
initially by a cascade lead filter with transfer function: 



G c (s) 



s + 2 
s + 20 



( 6 - 68 ) 
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The compensated system is shown in figure (6-27) . 



The characteristic equation of this system is: 



s 4 +26s 3 +125s 2 +(100+k) s+2k = 0 



(6-69) 



where k = 20. Equation (6-69) can be written in the form: 



s(s + l) (s + 5) C's+20) 

The Root locus of this system with k used as the 
Root locus variable parameter is shown in figure (6-30). 

At z = 0.6 the corresponding Root locus point is 
(£ = 0.60, = 1.232586). The value of w n although it 

could be measured directly from the Root locus diagram of 
figure (6-30) was found in the indicated accuracy by the 
use of the "singular point program". Therefore this 
initially compensated plant provides the desired singular 
root (assuming that the small difference between the desired 
value of w , i.e. , to =1.2 rad/sec and the actual value 
of = 1.232586 is within the system's tolerance). The 
shape also of the Root locus of figure (6-30) is such that 
promises some range of dominancy of the singular root. 

Based on these remarks the initially compensated 
plant was next compensated by a singular cascade compen- 



sator. the transfer function of which according to equation 
(6-18) is: 



k 



s + 2 



1 



(6-70) 




g- (s + 2) 



(6-71) 



s (s+1) (s+5) (s+20) +a (s+2) 
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The compensated singular system is shown in 
figure (6 -2 8a) . The characteristic equation of this 
system was found to be: 

s 8 +52s 7 + 9 26s 6 +C6700+a)s 5 +(.2Q825 + 2 8a)s 4 + 
(25000+177a)s 3 +(I0000+350a +208)s 2 + (6-72) 

(200a+808)s+808 = 0 

When the singular root (£ g = 0.60, w ns = 1.232586) 
was considered, the equation of the corresponding singular 
line was found from equations (.2-5) and (2-16) to be: 

8 = 3.31245248a - 219.4468287 (6-73) 

Introducing equation (6-73) into equation (6-72) 
and after some manipulations yields: 

s s + 28s 4 + 177s 3 +'416 . 2 49's 2 + 46 4 . 996s + 264 . 996 
a • . 

s 8 + 52s 7 + 926s 6 + 6700s 5 + 20825s tt + 25000s 3 + 5611. 0635s 2 

=-l (6-74) 

-17555. 746s- 17555.746 

or after factoring numerator and denominator yields: 

a : (s 2 + l . 479s + l . 519) (s + 20 .204) (s + 4.316> (s+2) 

(s 2 + l . 479s + l . 519) (s+20. 204) (s + 4.316) (s+19.786) 



(s+5.642) (s+1.4129) (s-0.840) 



(6-75) 
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which is in agreement with the conclusions made in para- 
graph four, i.e., the number of the fixed roots is four 
which is also the order of the characteristic equation of 
the system before the singular cascade compensator was 
introduced and that the Root locus of the compensated 
singular system with the parameter a used as the Root locus 
variable parameter has identical zeros and the same number 
of poles as the uncompensated system, i.e., the system 
before the singular cascade compensator was introduced. 

The Root locus of the compensated singular sys- 
tem was drawn by using a computer program in which equation 
(6-74) was implemented. This is shown in figure (6-31) . 

In this figure the four fixed roots are shown with dark 
circles at the corresponding locations. 

a. Range of Dominancy of the Singular Root 

By the use of a DSL simulation program in 
which the system of figure (6-28b) was implemented, the 
range of dominancy of the singular root (£ = 0.60, 

w ns = 1.232586) was investigated. More than one hundred 
properly selected parameter pair values were used within 
the system's stable region. The conclusion from this 
analysis was that the singular root has a very wide range 
of dominancy. The detailed results which were interpreted 
on both, the s-plane and the a-8 parameter plane, have as 
follows : 
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(1) S-plane . The Root locus of the com- 



pensated singular system is shown in figure (6-31) . This 
Root locus has been drawn with the parameter a - used as the 
Root locus variable parameter. According to equation 
(6-72) and (6-75), the compensated singular system has 
eight characteristic roots, four of which are fixed at the 
points : 

s 1 = -4.316 s 3 = -0 . 7 396+ j 0 . 9861 

s 2 = -20.204 s 4 = -0. 7396-jO. 9861 

i.e., two fixed real roots and the fixed complex pair of 
the singular roots. All these fixed roots are shown in the 
corresponding locations of figure (6-31) with dark circle 
marks. The location of each of the other four character- 
istic roots is varying on the corresponding Root locus 
segments as a function of the parameter a value. It must 
be noticed that a specific value of the parameter a defines 

a unique pair of parameter a and 3 values since for each 

Value of a corresponds one and only one value of 3 as can 
be seen from equation (6-73), which is the equation of the 
operating singular line. 

Interpreting the computer results which were 
obtained from a DSL simulation and a polynomial root finding 
program, the following remarks and conclusions can be 
stated in correlation of the parameter a value and the 
Root locus of figure (6-31) . 
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(a) 0 < a <66.25. Within this range of 

the parameter a value , the compensated singular system is 
unstable because there is a real root on the right half of 
the s-plane. Based on the need for a specific name to be 
given to this root, since it was found that after the 
singular root this is the only root which affects the 
transient response of the compensated singular system and 
the fact that this root lies on the same Root locus segment 
with the "singular root", it was called the "non singular 
root" . 



(b) a = 66.25. For this value of the 
parameter a the non singular root lies on the origin of 
tire s-plane and the compens-ated singular system is still 
unstable or conditionally stable. 

(c) 66.25 < a < 81.30. Within this range 
of the parameter a value, the non singular root lies on 
the real axis in the left half of the s-plane. The compen- 
sated singular system is stable and the non singular root 

is dominant. In figure (6-32) the , transient response of the 

■ t 

compensated singular system is shown for a unit step input 
and with an arbitrarily selected value of the parameter a 
within this range, i.e., a = 75 for which the corresponding 
value of 8 was found from the equation of the operating 
singular line, (equation C6--73)) to be 8 = 28.987. 
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(d) a * 81.30: For this value of the 

parameter a the non singular root is repeated at the 
point where the Root locus breaks away from the axis 
of reals, i.e., at this point the non singular root 
meets the other real root which was moving to the right 
on the real axis. This root is again the dominant root. 

In figure (6-33) the transient response of the compensa- 
ted singular system is shown for this value of the 
parameter a. As the value of the parameter a further 
increases, these two roots form a complex pair of roots. 

In this case by the name "non singular root" this complex 
pair will be meant. 

(e) 81.30 ^ a 95. Within this range 
of the parameter a value the non singular root is 
complex and although it has moved away from the real 

axis it is still closer to the origin and to the imaginary 
axis than the singular root. Within this range the con- 
tribution of the singular root to the transient response 
of the compensated singular system was negligibly small, 
and the non singular root was dominant. In figure (6-34) 
the transient response of the system is shown for a = 90. 
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(f) 95 v < a < 105. Within this range of 

the parameter a value the non singular root has ap- 
proached closer to the singular root. The contribution 
of both of these roots to the transient response of the 
compensated singular system is significant and none of 
these can be said to be dominant. 

When the second order model cor- 
responding to the singular root (? s = 0-6,o) ns = 1.2326) 

was considered, with transfer function: 



0) 



ns 



2 



2 

s +2 



r a) 

ns 



2 



S + 0) 



ns 



1.5193 

s 2 +0.7396s+1.5193 



(6-76) 



the corresponding transient response characteristics were 
calculated and they were found to be as follows: 



Time to peak of first 


overshoot : 


t = 3,19 
P 


sec 


Settling time 




t =5.41 
s 


sec 


Percent overshoot 




Mp =9.48 


% 


Transient oscillatory 


frequence : 


o)^ = 0.99 


rad/sec 
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Within this range of parameter a value the transient 
response characteristics of the compensated singular system 
were found to be as follows: 



Time to peak of first overshoot: 
Settling time : 
Percent overshoot : 
Transient oscillatory frequency: 



t = 5 sec. 

P 

5 «: t s <: 8 sec. 

8 £ Mp <: 111 

0.6 $: $ 1 rad/sec 



Because in practical applications the 
transient response characteristics which are of the basic 
interest are the settling time and the percent overshoot 
it seem from the above results that an acceptable solution 
to the problem can be derived within this range of the 
parameter a. In figure (6-35) the transient response of 
the system is shown for a = 100. 

(g) 105 < a < 131.5. Within this range of 

the parameter a value the singular and the non singular 
root are very close to each other. The contribution of 
both of these roots to the transient response of the com- 
pensated singular system is again significant, and none 
of these can be called dominant. The transient response 
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was in general characterized by long settling time, which 
was approximately tg =8.5 sec and high percent overshoot 
values, which were found to lie in the range from 11% up 
to 18.5%. In figure (6-36) the transient response of the 
compensated singular system is shown for a = 120. 



(h) a = 131.5. For this value of the 
parameter a the non singular root becomes identical with 
the singular root, i.e., the singular root is repeated. 
In figure (6-37) the transient response of the system is 
shown for this value of the parameter a. 



(i) 131.5 < a < 285. Within this range 

of the parameter a value the non singular root, although 
it is further away from the origin and the imaginary axis 
than the singular root, its contribution to the transient 
response of the compensated singular system is still 
significant and none of these roots can be called dominant. 
The transient response characteristics of the compensated 
singular system were found to be as follows: 



Time to peak of first overshoot 
Settling time 
Percent overshoot 
Transient oscillatory frequency 



2 . 8 < < 4 sec 

7<tg < 8.8 sec 
18.96 < Mp < 23.6% 

0.98 < < 1.1 rad/sec 



115 



In figure (6-38) the transient response of the compensated 
singular system is shown for a = 200. 



(j) 285 < a < 385. Within this range of 

the parameter a value the non singular root has moved 
further from the imaginary axis and from the origin, and 
the singular root has become the dominant root but with 
low degree of dominancy. Within this range, the transient 
response characteristics of the compensated singular system 
were found to be as follows: 



Time to peak of first overshoot 
Settling time 
Percent overshoot 
Transient oscillatory frequency 



6.6 < t < 7 sec 
" s - 

14.22 < Mp < 18.96% 

0.99 < £o t <1.1 rad/sec 



i.e., the percent overshoot is from 50% up to 100% over the 
corresponding to the singular root second order model 
value, and the settling time from 22% up to 30%. In figure 
(6-39) the transient response of the compensated singular 
system is shown for a = 330. 

(k) 385 < a < 1520. Within this range 

of the parameter a value the non singular root is so far 
away from the origin in comparison with the singular root, 
that its contribution to the transient response of the 
compensated singular system is either very little or 
neglegible and the singular root is the dominant root. 
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Within this range of the parameter a value the transient 
response characteristics of the compensated singular sys- 
tem were found to be as follows: 



Time to peak of first overshoot 
Settling time 
Percent overshoot 



11.85 $ Mp <? 14.221 
0.98 < 0 )^ <0.99 rad/sec 



2.8 < t <3.1 sec 



5.6 < t c < 6.6 sec 

v O s 



Transient oscillatory frequency 



i.e., the percent overshoot from 10% up to 50% over the 
second order model value and the settling time no more 
than 22% over that corresponding to the second order model 
value. In figure (6-40) through (6-44) the transient 
response of the compensated singular system is shown for 
several values of the parameter a within the range which 
was considered. It was noticed in this case that although 
the parameter values vary within a wide range, the transient 
response of the compensated singular system remains almost 
unaffected. 



of the parameter a value the non singular root appears again 
to affect appreciably the transient response of the com- 
pensated singular system due to the fact that although it 
is far away from the origin it is very close to the imagin- 
ary axis. The transient response is characterized with 
long settling time and high frequency of oscillations. In 
figure (6-45) through (6-47) the transient response of the 



(1) 1520 < a < 1922. Within this range 
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compensated singular system is shown for three selected 
values of the parameter a within the range which was 
considered. 



(m) a = 1922. For this value of the 
parameter a the non singular root lies again on the 
imaginary axis and the system is unstable or conditionally 
stable. The corresponding value of w for the non singular 
root at that point is w = 8.67 rad/sec. 

(n) a > 1922. Within this range of the 

parameter a value the non singular root lies on the right 
half of the s -plane and the compensated singular system 
is unstable. In figure (6-48) the transient response of 
the system is shown for such a value of a,i.e., a=2000. 

The conclusion from the above analysis is 
that there is in general a wide range of the parameter a 
values for which the singular root is dominant, and which 
may be restricted according to the tolerance of the 
transient response characteristic values which are usually 
determined with problem specifications. 

(2) Parameter Plane r In figure (6-49) the 
parameter plane diagram of the compensated singular system 
is shown for the singular values of £ s and u) ns together 
with the singular line which is associated with them. As 
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can be seen in this figure all the above mentioned lines 
are tangent at the same point of the parameter plane. 

In figure (6-50) the same parameter plane 
diagram is shown for different scale values. On the same 
diagram the compensated singular system stable region and 
dominant singular root region have been interpreted in 
terms of singular line segments. 

b. Hardware implementation of the compensated 

singular system 

According to Part C of this Section, the com- 
pensated singular system can be implemented as shown in 
figure (6-28a) or (6-28b) or (6-28c) , all of which are 
equivalent . 

Choice of the configuration to be used in 
practice will be in general determined from the specific 
application and the type of compensator available, i.e., 
integrated form or distinct components . 

When it is required to maintain the two para- 
meters a and 8 adjustable or when the compensator of the 
system is going to be built from distinct components then 
the configuration of figure (6-28b) seems to be the most 
representative one for such a case. 

Considering also the advance of the technology 
today in integrated circuits, it must be realized that 
any required compensator can be built in integrated form. 
In such a case either figure ^6-28a) or figure (6-28c) 
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can be implemented, after the two parameter values have 
been determined. 

c . Remarks 

When the singular line theory is considered as 
a design method for compensation of linear systems then 
the following remarks can be stated concerning the singular 
case which was considered in this section. 

Cl) It can be applied to any type of plant. 

(2) The design problem has been simplified and 
better organized since in general the basic design problem 
is not to find the required compensator but the appropriate 
values of the two parameters a and 3 which have to be used 
with the singular compensator which was defined by equation 
(6-18) . 

(3) The existence of the singular line pro- 
vides increased flexibility in the design procedure. 

It must also be noticed that when the two 
paramaters a and 8 are considered as time variant quantities 
of an adaptive system and its compensator respectively, then 
by adjusting their initial value and their rate of change 
in such a way that the operating point moves along the 
corresponding to the singular root dominant root region, 
system performance will be essentially unchanged. In 
this case other characteristic roots of the system will 
change but since they are not dominant, system performance 
can not be affected. When the above mentioned rate 
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adjustments are made automatically the compensated singular 
system becomes self-adaptive. This type of application 
of the singular line theory in adaptive systems needs to 
be further investigated. 

In the previous example the results 
which were obtained were related with a mathematical 
accuracy which in general can not be achieved in practice. 
In this case a reasonable question will be: How is the 

performance of the compensated singular system affected 
when the values of the two parameters a and 3 are set 
approximated with respect to their theoretical values? 
or what is the sensitivity of the singular dominant root 
with respect to each or both of the two parameter values? 
This problem is next analyzed. 

G. SENSITIVITY ANALYSIS 

1 . Introduction 

The characteristics of a physical component are 
subject to change for various reasons. Perfect control 
systems with stable components which can be built or can 
maintain their design values under any conditions are 
difficult to obtain in practice. Take for example such 
a passive component as a resistor. It is almost impossible 
to obtain a resistor with the exact design resistance 
value. For this reason a tolerance usually ranging from 
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1% up to 30% will accompany its specifications. Even if 
this "ideal" component is obtainable, factors like tem- 
perature variations, radiation etc., may eventually cause 
the resistance value to vary and produce undesirable 
results on the control system performance. These general 
rules are also applied when an active or solid state de- 
vice is considered. 

A measure of the effects of parameter variation 
upon the performance of a control system has in general 
been demoted as sensitivity. Basically there are two 
types of parameter variations studied by various inves- 
tigators. The first one which is the most popular and 
the most common in practice is the incremental variation 
case, where the change in parameter values is assumed to 
be very small, i.e., infinite-simal . The corresponding 
sensitivity study in this case is called <<Microscopic 
Sensitivity Analysis >>. The second case is that of 
relative large variations in parameter values i.e., 
finite changes. In this case the corresponding sensitivity 
study is called << Microscopic Sensitivity Analysis >>. 

The definition also of the sensitivity coefficient 
varies among investigators. Generally the sensitivity 
coefficient S can be defined as the ratio of the amount by 
which a control system’s performance characteristic deviates 
from its nominal (or original) value when one or more sys- 
tem parameters vary to its nominal value. From this 
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definition the sensitivity coefficient is mathematically 
given by: 



s = LS. 



( 6 - 77 ) 



where C denotes the nominal value of the system charac- 
teristic which is studied and Ac the deviation from this 
value when one or more system parameters which are con- 
sidered vary. 

There have been many sensitivity functions 
defined for the effects of parameter variations on the 
transfer functions of linear control systems, some of 
which are the following: 

a. The Kokotovic Method 

The root sensitivity functions developed by 
Kokotovic and Siljak provide most general solution to the 
sensitivity problem for small parameter variations in 
linear control systems and the most applicable to parameter 
plane methods. 

b. The F.H. Hollister Method 

The macroscopic root sensitivity functions 
defined by F.H. Hollister are applicable for large para- 
meter variations. 

c. Bode Method 

Bode defines a logarithmic or normalized 
sensitivity function of a variable, (t, d^), with 
respect to a parameter a^ , as: 
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s xi = 3 

aj 3 £n 

In general an analytic investigation of the 
sensitivity problem of a linear control system is neither 
simple nor easy. By computing sensitivity coefficients 
limits of individual parameter variations may be deduced 
and tolerance levels can be formulated. 

2 . Sensitivity Analysis of a Singular System 

The sensitivity analysis of a singular system 
can be performed by any method applicable to a linear 
feedback control system. By such a method the sensitivity 
of any system performance characteristic due to variations 
of one or more parameter values can be analyzed. Especially 
when compensation of a control system is designed by the 
use of the "singular line theory", then the transient 
response characteristics and the two variable parameters a 
and 3 are of the main interest for the singular system 
which results, at least for the first steps of the design 
process. Therefore it will be very important to know how 
transient response characteristics will be affected and 
how tolerance levels can be established concerning varia- 
tions of the two parameter a and 3 values. In terms of 
sensitivity this could be stated as: How- sensitive are 

the transient response characteristics of a singular system 
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to variations of the two parameter values and how can the 
corresponding tolerance levels be determined? 

The goal of the present analysis is to give an 
answer to these questions by a relatively simple and 
accurate way, in which computer use will be optional. 

Generally when designing by the use of 
singular line theory, specifications on transient response 
characteristics determine the desired operating singular 
point (see paragraph VI, F, 6). From this singular point 
or singular root the equation of the corresponding singular 
line will be evaluated. Some other specifications or 
physical limitations of the system may also restrict the 
permissible operating region on the singular line to a 
finite line segment. Then when the operating point moves 
or stays fixed on this singular line segment, performance 
requirements on transient response characteristics are 
satisfied. The design process up to this point has been 
described in detail in the previous part of this section. 

It will be assumed now, that the operating point 
moves away from the singular line due to variations of 
either one or both of a and 3 parameter values from their 
design values. This will be in general the case for an 
actual system. In such a case the acceptable tolerance 
on the transient response characteristics will determine 
how far away from the singular line the operating point 
is permitted to move. 
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Since there is at least one singular line for 
each value of £ for a singular system it will be expected 
that for a singular system there is an infinite number of 
singular lines on the parameter plane, each of which 
corresponds to a different value of C • Therefore accord- 
ing to the tolerances, which are given for the transient 
response characteristics of the compensated singular 
system, a permissible region can be established on the 
parameter plane for the operating point, which will be 
bounded between two such singular lines. This concept 
not only was applied in the following example but was 
also further investigated from a theoretical point of view 
and it was found that it is true and valid, which will be 
shown later in the example. 

Determination of these two singular lines on the 
parameter plane will not only give the answer to the pre- 
vious question, but it will give also a lot of other 
information concerning sensitivity and performance of a 
singular system for parameter variations. It will also 
give a precise way for tolerance values establishment for 

* t 

a singular control system. 

Although different specifications or different 
systems may give rise to different boundaries for the 
operating point on the parameter plane or different cri- ' 
teria in order to determine these boundaries, the procedure 
which must be followed will be in general the same. This 
will be illustrated by the following example. 
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3 . Example (6-7) 

Consider again the plant of example (6-6) with 
transfer function: 



G(s) 



20 

s (s + 1) (s + 5) 



which is compensated by an initial cascade compensator 
with transfer function: 



G c (s) 



s + 2 
s + 20 



The Root locus of the resulting system which is 
drawn with the gain of the plant used, as the Root locus 
variable parameter, is shown in figure (6-51). 

It is already known (see paragraph VI, F, 6) 
that any point of the above Root locus, which does not 
lie on the real axis is a singular point for the compen- 
sated singular system. It will be assumed that the 
specifications on the transient response characteristics 
for the compensated singular system require the operating 
point to be the Root locus point which is located at 
(? s = 0.60, w ns = 1.23) and that according to the given 
tolerance on the transient response characteristics which 
are interpreted in terms of the second order model, this 
singular point is permitted to move within the finite Root 
locus segment which is bounded between the Root locus 
points :. 
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0.65, w' 

’ ns 



= 1.03) 



c c; 



and 



( c" =0.55, w" = 1.58) 
v s ’ns 

These points can be determined by the intersection 

of the Root locus with the corresponding constant t, lines 

(see figure (6-51)) or can be evaluated by using the 

"Singular point program" in the computer. In the second 

case, although more precise values can be obtained for 

to ' and to " it will be seen later in this example that 
ns ns r 

such an accuracy is not actually needed. 

Since the above defined Root locus segment has 
been established by interpreting problem specifications 
in terms of second order model, will be actually the 
permissible region for the operating singular root, if 
and only if, this singular root has a high degree of 
dominancy anywhere within this region. 

By, either studying the Root locus of figure 
(6-51) or taking into account the results obtained from 
example (6-6) , a basic remark concerning the dominancy of 
the singular root in the case which is considered is that 
as decreases, the degree of dominancy of the singular 

root also decreases because it moves away from the origin 
and from the imaginary axis, and the contribution of the 
non-singular root (see example (6-6)) becomes more 
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significant. Based on this remark and the need for the 
given specifications on the transient response charac- 
teristics to be met, it was decided the permissible 
region for the operating singular root locus to be restric- 
ted between the points: 

(5=0.60, oi =1. 23) 
v s ns 



and 



(?' = 0.65, w' = 1.03) 

s ’ ns v 

which are shown in figure (6-51). 

At this point it must be noticed that such a 
decision is based absolutely on personal judgments the 
correctness of which can be checked by simulation of the 
compensated singular system at the limiting acceptable 
values for the two parameters a and 3, which were obtained 
from this analysis, as it will be shown later in this 
example . 

By introducing a singular cascade compensator, 
the characteristic equation of the compensated singular 
system becomes: 

s 8 +52s 7 +926s 6 +(6700+a)s 5 +(20825+28a) s 4 + 

( 2 5000+1 77a) s 3 + (10000+ 350a+ 208) s 2 + (6-78) 

(200a+808)s + 808 = 0 
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When the singular point (c = 0.60,w ns = 1.23) 
was considered, the equation of the corresponding singular 
line was found from equation (2-5) and (2-16) to be: 

3 = 3.304a - 218.355 (6-79) 

Comparison of. equations (6-79) and (6-73) reveals 
that the use of the computer for evaluation of the polar 
coordinates of a singular point and from these the equation 
of the corresponding singular line to be calculated, does 
not give any significant improvement on the accuracy of 
the equation of the singular line, for analysis and design 
purposes . 

When the singular root ( 5 ^ = o.65, = 1.03) 

was considered the equation of the corresponding singular 
line was found from equation (2-5) and (2-16) to be: 

8 = 2.444a - 138.097 (6-80) 

The above two singular lines were drawn on the 
parameter plane diagram which is shown in figure (6-52). 
After taking into account the results of example (6-6), the 
characteristic equation of the compensated singular system, 
and equations (6-79) and (6-80) of the boundary singular 
lines, the stable region A B C A b’c'a, and the region 
D E F G in which the varying singular root has a high 
degree of dominancy were approximately established on the 
parameter plane diagram of figure (6-52). 
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When the operating point lies at the point A 
where the two singular lines intersect, then the compen- 
sated singular system has both the corresponding singular 
roots as characteristic roots. When the operating point 
lies within the region D E F G where the varying singular 
root is dominant with a high degree of dominancy, then 
the specifications on the transient response characteris- 
tics for the compensated singular system are met, i.e: 

Peak overshoot Mp: From 6.8% up to 12.6% 

Settling time 't : From 4.6 sec up to 5.9 sec 

Trans, oscil. freq. : From 0.78 rad/sec up to 1.32rad/sec 

This was actually verified in the computer by a 
DSL simulation program of the compensated singular system, 
in which several pairs of a and 8 parameter values which 
define points in the region D E F G of the parameter plane 
diagram of figure (6-52) were considered. In figure (6-53) 
through (6-56) the time response of the system for some 
of these points is shown. 

a. Sensitivity Study 

Each operating point within the region D E F G 
of figure (6-52) lies on a singular line (it will be shown 
later in this example) , which corresponds to a singular 
root (? s ,to ns ). These polar coordinates of this singular 
root can be approximated for any such point from the 
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corresponding values of the boundary singular lines. For 
example the point and any other point which lies half 
way between the boundary singular lines will correspond 
to a singular root (C = 0.625, w ns = 1.13). When such a 
point, i.e., the point (a = 900, 3 = 2400) is 
considered the corresponding transient response charac- 
teristics can be evaluated from the second order model 
which corresponds to the singular root (C g = 0.625, w ns = 
1.13). 

Considering now that due to variations of 
either one or both of the parameter a and 3 values the 
operating point moves to another location within the 
region D E F G of the parameter plane of figure (6-52), 
then the corresponding new ? s and u) ns values can be 
approximated and from these the corresponding transient 
response characteristics. Therefore the sensitivity 
coefficient which relates the transient response 

i 

characteristics and the variations of the two parameter 

a and 3 values can be evaluated and from these sensitivity 

studies can be made as it is shown below: 

(1) Sensitivity of the Peak Overshoot due to 

Variations of the Parameter a . Let the 

operating point move from M^(a = 900, 3 = 2400) to M 2 

(a = 990, 3= 2400) due to a +101 change of the parameter a 

value. The corresponding singular points were approximated 

as (c = 0.625, w = 1.13) and (C ' = 0.64, id' = 1.07) 

^ s ’ ns s ns 
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respectively. The corresponding percent peak overshoots 
as they can be found from the second order model are 

r 

Mp = 81 and Mp =7.3%. Therefore according to the defi- 
nition of the sensitivity coefficient which was given in 
paragraph one: 



Mp = Mp* - Mp 
a Mp 



7.3-8 

8 



0.0875 



(2) Sensitivity of the Peak Overshoot due 
to Variations of the Parameter g . Let 
now the operating point more from M^( a = 900, 8 = 2400) 
to Mj(a = 900, 8 = 2640) due to a +10% change of the para- 
meter 8 value. The corresponding singular points were 

approximated as (£ = 0.625, w = 1.13) and (5 " = 0.61, 

s ns s 

= 1.19) respectively. The corresponding percent 

overshoot was evaluated from the second order model to be 
Mp = 8 % and Mp" = 8.9%. Therefore the sensitivity co- 
efficient which relates the peak overshoot and the + 10 % 
variation of the parameter 8 value from its value at the 
initial point M^ will be: 



Mp _ Mp - Mp _ 8.9-8 
’8 Mp 8 



0 . 1125 



Comparison of the above evaluated two sensi- 
tivity coefficients reveals that when the operating point 
M^ is considered, then the peak overshoot of the compensated 
singular system is more sensitive to variations of the 
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parameter 3 than it is to variations of the parameter a 
(almost 1.28 times more sensitive), and that the effect 
of the variations of these two parameters on the Peak over- 
shoot of the compensated singular system is opposite. 

It must be noticed that although the values 
obtained for the sensitivity coefficients are approximate, 
the conclusions which can be derived from a relative com- 
parison must be considered as very accurate for analysis 
and design purposes. 

In the same way any sensitivity problem which 
relates the transient response characteristics of the com- 
pensated singular system and variations of either one or 
both of the parameter a and g values can be studied. 

b. Selection of the Operating Point 

In general the operating point M of a compen- 
sated singular system must be selected where the singular 
root has the highest degree of dominancy (which results in 
the best possible approximation of the transient response 
characteristics of the compensated singular system by the 
corresponding second order model) and on a point halfway 
between the boundary singular lines (which will give the 
maximum tolerance for variations of the two parameters a 
and g, when it is considered that both can vary simul- 
taneously) . Under these considerations the best location 
for the operating point in the case which was studied in 
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this example seems to be at the point M(a = 1310 , 3 = 
3500), which is shown in figure (6-52) , It must be noticed 
that some other specifications of system's performance, or 
other limitations may require another location for the 
operating point M to be determined within the region 
D E F G. 



c. Tolerance Levels Establishment on Parameter a 

and 3 Variations 

After the operating point M has been chosen, 
then tolerance levels on the two parameters a and 3 
variations must be established. These tolerance levels 
must be determined in such a way that the operating point 
will be restricted to lie within a circle which has center 
the point M itself and is tangent to the two boundary 
singular lines and to any other boundary line in the 
neighborhood of this point. For the case of the example 
which was studied the required tolerance levels for both 
the a and 3 parameters are ±7.5% from the corresponding 
nominal values of a = 1310 and 3 = 3500, as it can be 
evaluated from figure (6-52), i.e., for no more than 7.5% 
variations towards both directions of each parameter a 
and 3 value, the operating point will lie within the circle 
shown in figure (6-52) whose center is at the point 
and the performance of the compensated singular system will 
be such that the specifications will be met. 
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4. Comments 



In establishing tolerance levels in practice some 
trade offs can be made by. considering each parameter 
separately since the nature of each is different, i.e., the 
parameter 8 is always associated with a gain value while 
the parameter a is associated either with another gain 
value or with pole locations (see figure (6-28a,b ,c) ) . If 
for example, in the previous case, which was studied above, 
the poles of the singular compensator can be set very 
accurately, due to the nature of the singular compensator 
which is going to be used, then its gain 8 will not be 
required to be set with a higher than ±15% accuracy, in 
order for the operating point to be restricted within the 
desired circle which was mentioned above, as can be seen 
from figure (6-52). 

Another general remark which was derived from 
figure (6-52) is that as far away from the intersection of 
the two boundary singular lines the operating point M 
is located, the less sensitive the compensated singular 
system becomes to the parameter a and 8 variations. 

The above comments together with the results 
obtained up to this point from the sensitivity analysis 
reveal that design of compensation by using singular line 
theory leads to very accurate setting of tolerance levels 
for the two parameter a and 8 values and very easy and 
accurate way of determining the optimal operating point for 
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the compensated singular system either on the s -plane or 
on the a- 6 plane. 

Although all the required information in order to 
make a conclusion about the sensitivity of a compensated 
singular system transient response characteristics, due to 
variations of the parameter a and $ values, have not been 
derived, based on the results of the analysis which was 
performed up to this point, it is believed that in general 
compensation of a system by using singular line theory 
may lead to control systems which are not sensitive or at 
least very sensitive to parameter values variations. 

In the previous example it was shown that when 
the operating point moves within the region D E F G of 
figure (6-52) then the transient response characteristics 
of the compensated singular system can be well approximated 
by the corresponding second order model. Actually there is 
no question about it, when the operating point lies on a 
singular line, like the two boundary singular lines D E 
and G F of the previous mentioned operating region. In the 
case where the operating point lies between these boundaries, 
although the results from example (6-7) reveal that the 
method is still valid, some reasonable questions which may 
rise, in this case from a theoretical point of view are: 

Is the compensated system still singular? Do singular 
conditions exist? Does the operating point still lie on 
a singular line? or have all the characteristics of the 
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singular conditions been lost? These questions were 
answered as follows: 

When a linear control system with negative unity 
feedback, whose open loop transfer function is GCs) = 



the characteristic equation of the compensated singular 
system will be: 



It has been shown (see paragraph III, C, 3) that such a 
system has in general an infinite number of singular points, 
i.e., points of the s-plane which satisfy equation: 



where n is an real integer or zero. When such a singular 
point has been specified, the equation of the corresponding 
singular line can be determined either by direct substitu- 
tion of the value of the singular point s in the character- 
istic equation of the compensated singular system or by 
using equations (2-5) and (2-16). From the equation of 
this singular line then, an infinite number of a and $ 
pair values corresponding to the singular point under 
consideration, can be determined. In the opposite case 
where a pair of a and B values has been specified the 




1 +a F (s) +B kF 2 (s) = 0 
g g 
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equation of the corresponding singular line and correspon- 
ding to this equation the singular point can also be 
determined . 

This will be illustrated by considering the follow- 
ing two cases from example (6-7) where the operating point 
moves within the region D E F G. In the first case the 
operating point arbitrarily was assumed to lie at the point 
M^(a=900, 6 = 2400) of figure (6-52). Then the characteristic 
equation (6-78) of the compensated singular system becomes: 

s 8 +52s 7 +926s 6 + 7 600s 5 +46025s 4 +184300s 3 + 373000s 2 
+372000S+192000 = 0 



By the use of a computer program with the ZPOLR 
subroutine the roots of the above equation were found to 
be : 



s 1 = - 0. 7044558 + j 0. 8843676 
s 3 = - 0.8651241+j5. 8923703 
s 5 = - 22.1206845 
s 7 = - 4.4148436 



s 2 = - 0.7044558-j0. 8843676 
s 4 = - 0.8651241-j5. 8923703 
s 6 = - 20.1762447 
s g = - 2.1490673 



Each of the above complex roots was checked and 
it was found that each lies on the Root locus of figure 
(6-51). Therefore each of these complex roots is also a 
singular root, (see statement after equation (6-61)). A 
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more general theoretical proof of this statement which 
is always true, can be derived from paragraph VI, F, 4. 

Because the operating region of the parameter 
plane, where the operating point lies has been designed 
as dominant root region only the first singular root, i.e., 
the singular root s i> 2 » which also was dominant, was further 
considered. This root is: 

S 1 2 = - 0.7044558 ±j 0.8843676 
or 

CC = 0.623, w „ = 1.131) 

v s ’ ns J 

When this singular root was considered, from 
equation (2-5) and (2-16), the equation of the corresponding 
singular line was found to be: 

3 = 2.8472a - 162.1855 (6-81) 

Since the coordinates of Mp i.e., a = 900 and 
8 = 2400 satisfy the above equation it implies that this 
point lies on a singular line of a dominant singular root 
whose coordinates are (C = 0.623, w = 1.131). 

Two points are interesting here. First, the above 
coordinates of the dominant singular root were approximated 
from the corresponding values of the two boundary singular 
lines of figure (6-52), during the sensitivity study of 
example (6-7), (see paragraph G3 a (1) of this section), 
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and second, when the other singular root was considered, 
then by the same method, it was found that the operating 
point lies also on the singular line corresponding to 
this root, i.e., finally the operating point lies on 
the intersection of two singular lines. Of course under 
dominant root considerations this last singular line is 
of no interest. 

It was next assumed that the operating point was 
moved to a new location Mg (a = 1430, g = 4280) of the 
dominant root region D E F G of figure (6-52). Then the 
characteristic equation of the compensated singular 
system was evaluated to be: 

s 8 +52s 7 +926s 6 +8130s 5 +60825s 4 +278110s 3 + 596100s 2 
+628400 s +'342400 = 0 

By the same method as in the previous case the roots 
of this equation were found to be: 



Sj^ = - 0.7257602 + j0.9470602 
s 3 - - 0.3981063 + j7. 5192493 
s 5 = -■ 23.1172585 
s ? = - 20.1934085 



s 2 = - 0.7257602-j0. 9470602 
s 4 = - 0.3981063-j 7. 5192493 
s 6 = - 4.3550711 
s g = - 2.0865288 



Each of the above complex roots was also checked 
and it was found that it lies on the Root locus of 
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figure (6-51). Therefore they are also singular roots. 

In this case the dominant singular root is: 

s ± 2 = - 0.7257602 ±j0. 9470602 
or 

(C = 0.608, w = 1.193) 

v s ’ ns J 

\ 

When this singular root was considered, from 
equation (2-5) and (2-16), the equation of the corresponding 
singular line was found to be: 

8 = 3.1307a - 196.1499 (6-82) 

Since the coordinates of the operating point Mg 
which was considered, i.e., a = 1430 and 8= 4280, satisfy 
the above equation, it implies that the operating point 
Mg lies on the singular line which is defined by equation 
(6-82) . It must also be noticed that equivalent remarks 
like these which were stated in the previous case about the 
other complex root s^ ^ also were found to apply and in 
this case. 

By the same way the singular lines and singular 
points corresponding to any point within the region D E F G 
or to any other point of the parameter plane, can be 
evaluated . 

The conclusion from these results is that indepen- 
dently of where the operating point lies on the a- 8 plane, 
singular conditions always exist. 
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The results from the example which was studied are 
summarized in Tabie I of Appendix G. 

H. STEADY STATE ACCURACY OF A SINGULAR SYSTEM 

1 . Introduction 

A measure of control system accuracy is usually 
the steady state error. Related concepts and definitions 
can be found in reference [2] . 

When steady state error considerations are made 
for a singular compensated system shown in figure (6-5), 
the question which arises is: What is the effect of the 

singular cascade compensator on the error constant K 
(which is also called the DC gain or zero-frequency gain) , 
of the system which is compensated by a singular cascade 
compensator? 

The goal of the present study is for an answer to 
be given to this question and for the required design 
procedure to be established. 

2 . Steady State Error Analysis of a Singular System 

It is known that if a finite error exists in steady 
state operation of a control system its magnitude is deter- 
mined by the reciprocal of the error constant K . 

A 
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When a unity feedback control system is considered 



with open loop transfer function: 



G C s ) = k F g (s) 



(s + z^ (s + z 2 ) . ... (s + z m ) 

k - | "■ - - - 

S (s+p i )(s+p 2 ) ... (s+p n ) . 



then the error constant K is: 

X 



K 



x 



£im s N G(s) 
s+o 



(6-83) 



or 



K = 
x 



k z z . . ,.z. 

1 2 m 

P P • . -P„ 
n 



(6-84) 



When the above considered system is compensated 
by a singular cascade compensator whose transfer function 
G sc ( s ) is given by equation (6-18), then the error constant 
of the compensated singular system is: 



K 



sx 



£im s N G (s) G(s) 
s->o sc 



(6-85) 



or 



K. 



sx 



kg z z . . . z 
i 2 m_ 

a p p ... p 
l 2 n 



K 

a x 



( 6 - 86 ) 
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i.e., the error constant of the compensated singular 
system differs from the error constant of the system 
before the singular cascade compensator has been introduced 
by a factor equal to the ratio 3/a of the two parameter 
values. Since an operating singular line does not in 
general pass through the origin of the a- 3 plane, the 
ratio 3/a will not in general remain constant when the 
operating point is moved along a singular line. 

When an operating region on the a-3 plane is 
considered, which actually will be the case, i.e., the 
parameter values will be fixed, then the error constant 
of the compensated singular system will have also a fixed 
value . 

Since the ratio 3/a corresponds to the slope of a 
straight line through the origin of the a-3 plane, it 
implies that the locus of all the points of the a-3 plane 
with a fixed value of the ratio 3/a is a straight line 
through the origin, and that the locus of all the points 
of the a-3 plane with a specified range of values of the 
ratio 3/a is a region of the a-3 plane which is bounded 
between two such straight lines through the origin. 

Therefore interpretation of any steady state 
accuracy specifications in terms of the 3/a ratio value 
leads to "steady state error boundary lines" establishment 
on the a-3 plane as is shown in the following example. 
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3 . Example (6-8) 

In this example it was assumed that steady state 
accuracy specifications were established for the compensated 
singular system of the previous example (6-7) , and that 
interpretation of these in terms of the error constant 
K set the following restriction on the 8/a ratio value, 

b A 

i . e . : 



> 1 (6-87) 

The corresponding "steady state error boundary 
line" for this case is the straight line through the origin 
with slope equals to one. This was drawn on the a- 8 plane 
shown in figure (6-52) . Then all the points of the a-8 
plane above this line satisfy steady state accuracy 
specifications of the problem. Therefore any point of the 
dominant singular root region D E F G shown in figure (6-52) 
will also give an acceptable steady state error constant 
for the compensated singular system. 

In the same way, when another system or different 
specifications are given concerning the steady state 
accuracy establishment of the corresponding "steady state 
error boundary lines" or line, on the parameter plane the 
required information for analysis and design purpose can 
be derived. 
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VII . 



THIRD SINGULAR CASE ANALYSIS 



In this section the "third singular case" which was 
defined in Section IV, was considered for analysis. This 
singular case corresponds to the compensation structure of 
figure (4-5) where: 



H(s) = as 

and the plant G(s) = kFg(s) has at least two poles at 
the origin, i.e., the cascade compensator has no zero's 
at the origin. 

The characteristic equation of the compensated 
singular system is: 

1 + aks F (s) + Bk (sF (s) 2 ) = 0 

S § 

A. SYSTEM OUTPUT AND ERROR ANALYSIS 

The transfer function of the compensated singular 
system is: 



G (s) = kF g (s) 
Gc( s ) = F g (s) 



(7-1) 



T (s ) = £IiL_ = 



eks 2 Fg 2 (s) 



(7-2) 



R(s) 
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and the output and error are respectively; 
C(s) = R(s) TCs) 



C7-3) 



and 

E(s) - R(s) {1 - T(s) } (7-4) 



1 . Plant Type Two 

When the plant is type two, i.e.: 



G(s) = kF g (s) = k- 



(s+Zj) (s+z 2 ) . . .(s+z m ) 
s 2 Cs+p 1 ) (s+p 2 ) . . . (s+p n ) 



(7-5) 



equations (7-3) and (7-4) become respectively: 

3k [~~ |(s + Zj) 2 



C(s) = 



3 = 1 



n m n m 

/n (s+p k ) +aks I |(s + Zi) I |(s+pj c ) + 3 k n (s+Zj ) 
k=i j = i k= l i = l 



-R(s) 



3 = 1 



(7-6) 



and 



E(s) = 



n m n 

f ICs+Pk)' + aks \ \(s + Zi ) | l(s+p k ) 

k= l j = i k= l 



m 

l ( 



-R(s) 



|(s+Pi c ) 2 +aks ( |(s + zj) i i(s+pk)+ 3 k } | (s + Zj) 

3=1 k= l j=i 



k= i 



(7-7) 



When the final value theorem was applied to the 
above equations the steady state output and error were 
evaluated assuming that all the poles of the transfer func- 
tion T(s) lie on the left half of the s-plane. 
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a. Step Input 



When a step input R(t) = AuCt) or R(s) » A/s 
was considered it was found that: 



C(t) = C ss Ct) = A 

t -HX ' 

and 

ECt) - E ss Ct) * 0 

t-** 



b. Ramp Input 

When a ramp input R(t) = BtuCt) or R(s) 

was considered it was found that: 

CCt) = C ft) =« 
t+“ 



and 



E(t) - E (t) 

t-J-CC 



Ba 



8 



n 

k= l 
m 

n 



p 



k 



z . 
3 



B/s 



2 



R(s) 



and 



c. Parabolic Input 

2 

When a parabolic input R(t) = Ct u(t) or 

3 

2C/s was considered it was found that: 

c(t) = c ss (.t) = - 

t->oc 

E(t) = E ss Ct) = cc 

t-Xa 
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2 . Plant Type Three 



In the same way, the steady state output and error 
value for a type three plant were found to be as follows: 



and 



and 



a. Step Input 

C(t) = C ft) - A 

t+ <x 



ECt) = E ss Ct) - 0 

b. Ramp Input 

CCt) - C Ct) - « 

t + « 



E(t) = E ss Ct) 
t-> « 




and 



c. Parabolic Input 
CCt) = c Ct) = oc 

t->oc 



E(t) = E ss (t) = ex 

t-*-<x 
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3. Remarks 



The general form of the characteristic equation 
of the compensated singular system in the case which was 
studied in this section is: 

1+ctks F,Cs) + l?k{s F Cs)} 2 = 0 (7-8) 

& o 

where the transfer function of the plant is: 

G (s) - kF g (s) (7-9) 

By defining: 

F (s) = sF ff (s) 

C> 

the characteristic equation (7-8) becomes: 

1+akF (s) + gk F 2 (s) = 0 (7-10) 

which is the characteristic equation of a singular system. 

Comparison of equations (7-10) and (6-16) reveals 
that the design process which has to be followed when 
compensation is performed by the compensation scheme, 
which corresponds to the "third singular case" will be in 
general the same with that corresponding to the "second 
singular case”. The only difference will be that instead 
of using the Root locus of the uncompensated system as the 
basic design tool, the Root locus of the uncompensated 
system with one pole less than the actual number of poles 
at the origin has to be used. This will be illustrated in 
the following example. 
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B. EXAMPLE (7-1) 



In this example the plant shown in figure (7-1) was 
considered with open loop transfer function: 

G (s) * kF C s ) « 100 -r^ (7-11) 

g s Z Cs + 2) 

This plant is unstable and it was assumed that prob- 
lem specifications dictate for a stable system with a 
dominant root with polar coordinates: 

C = 0.40 

and 

2 < w <3 

n 

The compensated singular system is shown in figure 
(7-2) . As can be seen the whole design problem consists 
of finding the appropriate pair of a and 3 parameter 
values by which the given specifications can be met. This 
was performed in three steps which follow: 

1 . Step One 

1 f 

The first step was to determine the required 
operating singular root for the compensated singular sys- 
tem. 

The Root locus of the uncompensated system with a 
pole less than the actual number of poles at the origin 
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is shown in figure (7-3), i.e., this is the Root locus of 
equation : 

k - *--l (7-12) 

s (s+2) 

By the same reasoning as in the case of the "second 
singular case" it can be shown that any point of the above 
Root locus, which does not lie on the real axis is a sin- 
gular point of the compensated singular system. (See Section 
VI, F, 4, b) . Considering next the specifications the 
desired singular point was located on the Root locus at the 
point (c = 0.40, w n = 2.5), 

2 . Step Two 

The second step was to find the equation of the 
singular line, corresponding to the selected singular 
point. 

From equation (7-9) the characteristic equation 
of the compensated singular system was found to be: 

s 4 + 4s 3 +(4 + 100a)s 2 + 200aS + 100j3 = 0 (7-13) 

When the singular root (c s s 0.40, w ns - 2.5) of the 
above characteristic equation was considered, the equation 
of the corresponding singular line was found from equation 
(2-5) and (2-16) to be: 



6 = 6. 25a-0. 390625 



(7-14) 



3 . Step Three 



The third step which was the last one was to find 
the required operating point of the singular line which 
was defined by equation (7-14) . 

Because it was also desired that the design pro- 
cess for this problem be completed without using a computer 
this step was performed as follows. 

Introducing equation (7-14) into (7-13) the 
characteristic equation of the compensated singular system 
becomes : 

s 4 + 4s 3 +(4 + 100cOs 2 + 200cts + 625ct-39.0625 = 0 (7-15) 

Although according to the results which were stated 
in Section VI, F, 4, the Root locus of the compensated 
singular system, i.e., the Root locus of equation (7-15), 
with the parameter a used as the Root locus variable para- 
meter has the shape of the Root locus of equation (7-12) 
shown in figure (7-3), and therefore the drawing of the 
Root locus of equation (7-15) is not actually required, it 
was done in this example for a more detailed illustration. 

Equation (7-15) can be written as: 

100s 2 + 200s + 625 _ 

a — - - 1 

s 4 +4s 3 +4s 2 -39.0625 



154 



or after factoring numerator and denominator; 

{ s+ (1+ j 2 . 291) } { s+ (1- j 2 ♦ 291) } ^ 

(s+(l+j2. 291) } (s+Cl-j2.291)} {s+3.693} {s-1.693} 

The Root locus of the above equation is shown in 
figure (7-4). This root locus which was drawn by hand, 
consists of a fixed point there where the singular root 
(S s =0.40, a) ns = 2-5) or s--l±j2.291 is located and the 
regular Root locus segments. 

Comparison of figures (7-3) and (7-4) reveals that 
the two Root loci have actually the same shape. 

Because the singular root and the complex non- 
singular root (See Section VI, F, 4 , a, (1) , (c) ) , have always 
equal real parts, in this case, in order for the singular 
root to be dominant its residue must be much greater than 
that corresponding to the non-singular root, i.e., the 
singular root must lie much closer to the origin of the 
s-plane than the non-singular root. Therefore the non- 
singular root must lie far away above the singular root 
which is shown on the Root locus either of figure (7-3) or 
(7-4). 

By arbitrarily selecting the distance of the non- 
singular root from the origin to be relative very large, 
in comparison with that which corresponds to the singular 
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root, dominance of the singular root in general shall be 
expected, e.g., let the non-singular root be: 

s = -1 + j50 

Then the left part of equation (7-15) must be equal to the 
product of the singular root and non-singular root. This 
yields a = 25.07. From equation (7-14) the corresponding 
value of the parameter 3, was found to be 3 = 156.31. 

Therefore the required parameter pair values or the 
desired operating point of the a-3 plane is (a=25.07, 
3=156.31) . 

4 . Remarks 

The compensated singular system shown in figure 
C7-2) was simulated in the computer by using a DSL program 
in which the two parameters a and 3 were set at the above 
indicated values. 

From this program it was found that the time 
response characteristics of the compensated singular system 
were exactly the same as these corresponding to the singular 
root second order model. It is shown in figure (7-5). 

It must also be noticed that the desired perform- 
ance of the compensated singular system can be achieved by 
an infinite number of a and 3 parameter pair values, which 
correspond to a different selection in the location of the 
non-singular root on the s-plane. 
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Criterion for the selection of the appropriate 
location for the non-singular root can be based either on 
sensitivity or and steady state accuracy considerations as 
follows : 

a. Sensitivity Considerations 

As the location of the non-singular root moves 
away from the axis of reals on the Root locus of figure (7-4) 
the value of the Root locus variable parameter a increases. 

Considering next the parameter a-3 plane and 
the singular line (equation (7-14)), which corresponds to 
the operating singular root Cs s = 0.40, w ns = 2.5), as the 
value of the parameter a increases the operating point 
(a, 3) moves away from the origin on the singular line, but 
according to the results obtained from the sensitivity 
study in the previous section it yields a less sensitive 
compensated singular system with respect to variation on 
the parameter a and 3 values. These remarks can be used 
accordingly in each case. 

b. Stealy State Accuracy Considerations 

When the error constant of the compensated 
singular system was considered it was found to be: 

K = JL. 

sx 

la 

When problem specifications concerning the 

error constant K are also established, then an additional 

sx 
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criterion for the selection of the location of the non- 
singular root will be the desired value of the parameter a 
for which the ratio 3/a is such that the optimal error 
constant coefficient can be achieved. 

In general the ratio 3/a can be effected either 
by selecting a different location for the non-singular root 
or by selecting a different operating singular point, i.e., 
operating on a different singular line or by introducing 
additional gain adjustments at the forward path of the 
minor loop of figure ( 7 - 2 ) . 

It must be noticed also that any other specifica- 
tions can be taken under consideration in each design step 
of the previous illustrated example. 
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VIII. 



EXTENSIONS OF THE SECOND SINGULAR 



CASE COMPENSATION STRUCTURE 
A. INTRODUCTION 

Among the three singular cases, which were previously 
analyzed, the most interesting for applications in design 
of compensation of a linear feedback system by the 
dominant mode method, appears to be the "second singular 
case" (See Section VI) . 

The basic advantages of this case are that the plant 
type number is preserved, it can be used for any type of 
plant and for any plant structure. 

On the other hand the basic limitation of the "first 
singular case" CSee Section V) is that in general the type 
number of the plant is not preserved, and for the "third 
singular case" (See Section VII) is that the plant has to 
be at least type two. 

For the above reasons, the "second singular case" was 
further considered and some extensions of its compensation 
structure were studied. 
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B. SECOND SINGULAR, CASE COMPENSATION STRUCTURES 



1 . The Basic Structure 

The basic compensation structure of the "second 
singular case" is shown in figure (6-4) where: 



G (s) = kFg(s) = k- 



(s + z x ) (S + Z 2 ) . . . (s + z m ) 

S (s + p^) + f S + P r( ) 



( 8 - 1 ) 



represents the transfer function of the plant to be com- 
pensated. 

H(s) = a (8-2) 

represents a gain adjustment, introduced as a negative 
feedback from the output to the input of the plant, and 



G c (s) = 3F g (s) 



( s + z 1 ) (s + z 2 ) . . . (s + z m ) 
s N (s + P 1 ) (s+p 2 ) . . . (s+p n ) 



(8-3) 



represents the transfer function of a cascade compensator 
whose poles and zeros are identical with these of the 
plant and its Root locus gain is the second adjustment of 
the compensation structure. 

This compensation structure defined above, cor- 
responds to a singular system which was studied in Section 
VI as the "second singular case". 
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Its transfer function has the general form: 




(8-4) 



The error constant 

was found to be: 

m 



K n 
sx-B 






K. r, was evaluated and 

SX “ D 



it 



(8-5) 



The characteristic equation also of this basic 
singular compensation structure was found to be in the 
general form: 

1+cdeF (s) + 6kF g 2 (s) = 0 (8-6) 

From this basic structure six modified singular 
compensation structures were derived. They were analyzed 
as follows: 
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2. Modified Structures 



a. First Modified Structure 

The first modified singular compensation 
structure is shown in figure (8-1) . It was derived from 
the previous basic structure and has been analyzed in 
Section VI. 

In figure C8-1) G(s) denotes the transfer 
function of the plant to be compensated and G sc (s) the 
transfer function of the singular cascade compensator, 
which according to equation (6-18) its transfer function 
is : 



3 



G sc Cs) - .— 1=1 



m 

( I (s+Zj 



) 



n m 

s^ 1 |(s + P ) +a 

k-i K 



n cs+z. 
j-i j 



(8-7) 



where s^, p^ and z - , represents the poles at the origin, 

the other poles and the zeros of the plant respectively. 

‘ ( 

The transfer function of this compensation 
structure was found to be: 



( 8 - 8 ) 



Tj(s) = 



C (s) 



m 

3k I l(s + z .) 
i-i J 



R(s) 



2N 



n n m m 

l(s+P ) 2+ as N ["^(s+P ) 1 j(s+z,)+3k | ](s+z. 

k= l K k=i K j=l J 3 = 1 J 
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The error constant K was also found (See 

sx-1 



equation (6-85)) to be: 



m 



K 




(8-9) 



sx-1 a 



a x 



n 




where K denotes the error constant of the plant, before 

A 

it was compensated by the singular cascade compensator. 

The characteristic equation of this compensa- 
tion structure will be in the general form: 



compensation structure with respect to the previous basic 
one, are the following two: 



structure is simpler because only a single cascade filter 
is required. Historically it is realized that series 
compensation has played a prominent role in system design. 

(2) Greater Error Constant . Comparison of 
equations (8-5) and (8-9), reveals that the error constant 
of this structure is greater than the error constant of 
the previous basic one by a factor equal to the Root locus 
gain k of the plant. 



1 + aFg (s)+6kF g 2 (s) = 0 



( 8 - 10 ) 



The basic advantages of this modified singular 



(1) Simpler Structure. The compensation 
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b. Second Modified Structure 

The second modified singular compensation 
structure is shown in figure (8-2). It was derived from 
the previous one by changing the location of the singular 
compensator from the forward path at the feedback path. 

Of course this modified singular structure is no longer 
a unity feedback control system. 

For this compensation structure the transfer 
function was evaluated and it was found to be: 



T 2 (s) = 



CCs) _ 
Ws) 



n m m 

ks N n (s+P ) n (s+z.)+ak j { (s+z.) 

k=i K j=i 3 j = i 3 

n . n m m 

s N | l(s+p ) +as N 1 |(s+p )| 1( s+z-)+3k i 1 (s+z.) 

k=i K k=i K j=i 3 j=i 3 



( 8 - 11 ) 



The characteristic equation also of this 
compensation singular structure was found to be: 

1+G(s) g'(s) = 0 



or 

l+<xFg-(s)+gk F g 2 (s) = 0 (8-12) 

In order for steady state accuracy to be in- 
vestigated the following considerations were made: 

The error of the system is: 



E=R-C=(1-T 2 )R 



(8-13) 
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Introducing equation (8-11) into equation (8-13) 



yields : 



E (s) 




m 

) + (3k-otk) n (s+z.) 
i-1 l 



2N 



n 

n 

k=i 



(s+p ) +as N 



n m 

n (s+p ) | 1 (s+z .)+3k 

k=l K j=i J 




2 



R(s) 



(8-14) 



When a type one plant was considered and a ramp 

2 

input, i.e., N=1 and R(s)=B/s , and when the final value 
theorem was applied to equation (8-14), assuming a stable 
system, the steady state error was evaluated as: 

E(t)=E ss (t) = - 
t-+« 

For a physical interpretation of the above re- 
sult the following considerations were also made on the 
compensation structure which is shown in figure (8-2). 

For a type one plant, and a step input, at 
steady state E has to be zero, i.e., E g = 0. Because 

E = R-C , it implies that R * C . Since C is a 
* r ss ss ss 

constant and the system is not a unity feedback system 

C ss * C ss- Therefore c ss * R - 

When a type one plant was again considered 

with a ramp input, at steady state E has to be constant, 
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i.e., E sS = Constant. Because E = R-C this implies that 
• • > 

R = C . Since C is a ramp and the system is not a 

3 3 S S 

• • > • • 

unity feedback system, C„ f C „ . Therefore C = R, 

which implies that an infinite error E gs must be expected. 
This result was also derived previously by analytical 
study of equation (8-14) . 

Based on these remarks and results, the con- 
clusion is that such a system does not reproduce the input 
signal at the output or in other words at steady state the 
output does not follow the input, which may be an advan- 
tage or a disadvantage depending on the specific applica- 
tion, and which can be eliminated as will be seen by the 
following modified singular compensation structure. 

When a type one plant was considered again and 
a step input applied, the steady state error was found to 
be : 

E ss= C 1 - f ) A 

where A is the value of the step input, i.e., at steady 
state the output differs from the input by the value of 

Ess- 

In order for all the above mentioned problems 
to be eliminated when this is desired, the following 
considerations were made on this modified singular 
structure . 
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When R was considered as the input signal which 
the output of the system was desired to reproduce, i.e., 
to follow, then the error will be given by equation: 

E = R-C (8-15) 

By considering also an input which is not followed 
by the output to be R , then according to the notations 
used: 

C = T 2 r" (8-16) 

Introducing equation (8-16) into equation (8-15) yields: 

E = R-T 2 r" (8-17) 

By assuming that the actual input to the system 

* 

R is correlated to the desired input R by a function Y, 
i.e: 

r"= YR (8-18) 

then equation (8-17) becomes: 

E = R-T 2 YR = (1-T 2 Y) R (8-19) 

Introducing equation (8-11) into equation (8-19) yields: 



E(s) 




(8-20) 
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In order for the steady state error to be 
finite for a ramp input to a type one plant or zero for a 
step input, Y has to be defined as: 

Y = £ (8-21) 

a 

Introducing equation (8-21) into equation (8-18) yields: 

r'-£r C8-22) 

a 

After the above result was obtained a third 
modified structure was derived and is analyzed as follows: 

c. Third Modified Structure 

The third modified singular compensation 
structure is shown in figure (8-3), where R is the input to 
the system which is to be followed at steady state by the 
output of the system. 

The transfer function of this modified struc- 
ture was found to be: 




(8-23) 
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The characteristic equation was also found to 



be in the general form: 

1+a Fg(s)+3k Fg^(s) = 0 



(8-24) 



By the same procedure as in the case of the 

first modified structure, the error constant K , for 

’ sx- 3 

this case was found to be: 

m 



K 

sx- 3 




a 3 

— r 

a - 3k 



K 



x 



(8-25) 



where K denotes the error constant of the plant before it 

was compensated. For the examples which were worked out 

2 

in Section VI, the ratio aS/{a -3k} turns to be greater 
than unity. Therefore, when such is the case the error 
constant of the uncompensated plant will be improved when 
the plant is compensated. 

When this compensation structure was compared 
with the first modified structure it was found that they 
have the same characteristic equation but different error 
constant and transfer function. Therefore the stable 
region on the a-3 plane will be the same in both cases. 
Also for the same operating point C a >$)> although they 
will have the same characteristic roots, the residues of 
these roots will be different in each case due to the dif- 
ferent numerator of their transfer functions. Since the 
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dominancy of a characteristic root is a function of its 
residue, in addition to the distance of the location of 
this root on the s-plane from the imaginary axis, a dif- 
ferent degree of dominancy for the operating singular 
root must be expected in general in each case. 

For a specific plant and a specific operating 
point on the a-3 plane a further investigation of the 
dominancy of the operating singular root can be performed 
either by a computer simulation of each compensation - 
structure and comparing their time response characteristics 
with those corresponding to the singular root second order 
model or by analytical evaluation of the residues cor- 
responding to the characteristic roots of each case'. 

d. Fourth Modified Structure 

The fourth modified singular structure is^ 
shown in figure (8-4). It was derived from the first one 
by properly defining the feedback compensator in order for 
these two compensation structures to have the same open 
loop transfer function, i.e.: 



G(s) 

GCs) G (s) = * (8-26) 

1+G(s) G (s) 
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from. which the transfer function of the filter G sc (.s) was 
evaluated to be: 



G sc< s > 



N 



n 



n 



m 



m 



I |(s+p k ){s N | l(s+p k )+a | l(s+z.)-B 1 ;(s+z.)} 

k=l K k=l K j-i 3 j-l 3 



m 2 

6 k j |(s+Zj) 

j = 1 



(8-27) 



Because G(s) defined by equation (8-1) denotes the 

transfer function of a physical control system, it implies 

that N+n>m, i.e., the above defined filter has more zeros 

than poles. In practice a structure of this type is 

realized by introducing additional poles as needed, but 

placing them well outside the bandwidth of the system. 

Whether this can be done without affecting the singularity 

characteristics has not been investigated. Based on the 

argument requirement for singular conditions (See equation 

(3-12 and (10-11)), it is believed that if these additional 

poles are placed on the s-plane in such locations that 

* 

their net contribution to the argument of G gc (s) is 
negligibly small, then it can be assumed that the singu- 
larity characteristics of the system remain unaffected. 
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The transfer function of this structure was found 



to be : 



t c s i - G(s) _ G(s) 

1 4 '- S J Rf s') -* 

KLSj 1+G(s)+G(s) G sc (s) 



or after some manipulations: 



(8-28) 



T 4 (s) 



- _CCs) - 
■W 



m 

8 k j | (s+z .) 2 

J zi 

n n m ’ m 

= 2N ri(s + p k ) 2 ^s N | ](s+p k ) I |(s+z.)+0kj |(s+z • ) 

k=i K k=i K j=i J j=i J 



(8-29) 



The characteristic equation of this compensation 
structure was found to be in the general form: 



l+aFg(s)+8k 



Fg 2 (s) = 0 



(8-30) 



By the same procedure also, as in the case of the 

first modified structure, the error constant K . for 

sx-4 

this case was found to be: 



K A = - K 
sx-4 a x 



(8-31) 



Comparison of the results which were obtained for 
the first modified compensation structure with these 
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obtained above, reveals that the performance of the system 
in both cases is the same. 

e. Fifth Modified Structure 

The fifth modified singular compensation 
structure is shown in figure (8-5). It was derived from 
the second one shown in figure (8-2) by splitting the 
feedback path into two parallel paths as shown, such that 
their sum is the compensator transfer function: 

G sc (s) 5 1+G s C (s) (8-32) 

from which the transfer function of the minor loop feed- 
back filter was found to be; 



m 



N 



n 

n 



,,, , CS-oO^'Cs+zp-s > ks \ (s*p k ) 
sc L J 



n 



m 



S N I |(s+p k )+a| | (s+Z.) 

k=i K j=i J 



(8-33) 



The transfer function of this compensation 
structure is: 



T 



s 



(s) 



C(s) = G(s) 

RCs) 1 + G(s)+G(s) g” sc (s) 
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or after some manipulations it becomes: 



T s (s) 



. C(s) . 

~Ws J 



n m m 

ksN | l( s+ P k ) [Hcs+zJ+okj" - ] (s+z,) 2 
k=i j=i J j = i J 




n m m 

n C s* Pv ) r j(s+z.)+gk i i(s+z.) 2 

k=i K j=i J j=i 3 



(8-34) 

The characteristic equation of this compensation 
structure was found to be in the general form: 

l+aF g (s) + gkF g 2 (s) = 0 (8-35) 

Also when steady state accuracy considerations 
were made for this compensation structure the results were 
the same as these for the second modified structure. Com- 
parison of the results corresponding to these two cases 
reveals that they are equivalent under performance con- 
siderations . 

f. Sixth Modified Structure 

This is shown in figure (8-6), and it was de- 
rived from the third modified structure with the same 
reasoning by which the fifth modified structure was de- 
rived from the second one. 
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It was also found that the third and sixth 
modified structures have the same transfer function, 
characteristic equation, error constant and thus under 
system performance considerations they are equivalent. 

3. Remarks 



Although the "second singular case" can be used 
for compensation of any control system, i.e., it may have 
general application, its "basic structure" has certain 
disadvantages, which have been stated in the previous 
analysis and which may not be acceptable in some cases. 

When the characteristic equation- of the basic 
structure is considered, it will have the general form: 

1+a.kF (s)+gkF 2 (s) = 0 
& S 

On the other hand the six modified structures, 
which were analyzed above, have certain advantages over 
the "basic structure" and a common characteristic, i.e., 
they have the same characteristic equation, which has the 
general form: 

1+cxF (s)+gkF 2 (s) = 0 

o o 

Also their performance characteristics can be mathematical- 
ly correlated. 

This fact aids the designer in making his decision 
on the problem of the structure without additional design 
effort and protects from arbitrary choices. 
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In the following example all these modified 
singular compensation structures were used to compensate 
the same plant in order that the derived analytical 
results could be compared. The compensated system per- 
formance was also studied in each of these cases. 

4 . Example (8-1) 

In this example the plant of example (6-7) was 
considered for which the transfer function G(s) of the 
system before any singular filter has been introduced is: 

20 (s + 2) 

G (s) = (8-36) 

s (s + 1 ) (s + 5) (s+20) 

Also as operating point from the dominant root 
region of the a - 3 plane corresponding to the singular root 
(t s =0.60, w =1.23), the point (a = 920, 3 = 2500) was ar- 
bitrarily selected. Then each of the six modified sin- 
gular compensation structures was considered separately 
and the performance of each for a step input, R=1 Volt and 
a ramp input R=0.1 rad/sec, was studied in the computer as 
follows : 



a. First Modified Structure 

This is shown in figure (8-1) . According 
to equations (8-7) and (8-36) the transfer function of 
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the singular cascade compensator G sc (s) was found 
to be : 

2500 (s+2) 

G sc (s) * 

s (s +1) fs + 5) (s+20) +920 fs+2) 



By simulation of this system in the computer 
with the use of a DSL program the time response character- 
istics were obtained for a step and a ramp input. They 
are shown in figure (8-1A) and (8-1B) respectively. 

The characteristics of the transient response 
corresponding to the singular root (? s = 0.60, cj ns = 1.23) 
were evaluated and they were found to be: 



t = 3.19 sec 
P 



t = 5.41 sec 
s 

Mp = 9.48% 

oj . = 0.99 rad/sec 



Time to peak of first overshoot 
Settling time 
Percent overshoot 
Transient oscillatory frequency 

When the step input was considered then it 
was found from the computer print output that: 

Steady state output value 
Time to peak of first overshoot 
Settling time 
Percent overshoot 
Transient oscillatory frequency 



C ss « R = 1 Volt 



tp 



= 3.2 sec 



t s = 5.4 sec 



Mp = 11% 

= 0.98 rad/sec 
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When the ramp input was considered then it 
was found also from the computer print output that: 

Steady state output value : C gs = R-0.002 

Steady state output rate : C gs = R = 0.1 

Settling time : ts =5.4 sec 

These results can also be derived approximately from the 
corresponding figures of the time response characteristics, 
figure (8 - 1A, B) . 

b. Second Modified Structure 



When the singular compensator G sc (s) was 
located at the feedback path as shown in figure (8-2), the 
transient response characteristics which were obtained are 
shown in figure (8-2A) and (8-2B) for a step and a ramp 
input respectively. 

The exact transient and other characteristic 
values which were read from the computer print output have 
as follows: 



(1) Step Input. 



Steady state output value 
Time to peak of first overshoot 
Settling time 
Percent overshoot from C 

ss 

Transient oscillatory frequency 



C = 0.368 = | R 
ss 3 

t =3.2 sec 

t = 5.4 sec 

Mp = 10.2% 
top = 0.92 rad/sec 
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(2) Ramp Input . 



Steady state output value : C ss followed a ramp 

different than the input ramp. 

• • 

Steady state output rate : C = 0.036 = ^-.R 

S S p 

Settling time : t =5.4 sec 

c. Third Modified Structure 

When a gain adjustment equals with 2.717, which 
corresponds to the ratio 3/a value, i.e: 

£ . _25£0 . 2.717 
“ 920 

was introduced at the indicated in figure C8-3) location 
the third modified singular structure was obtained. 

The transient response characteristics which 
were obtained in this case for the same step and ramp 
input as in the previous cases, are shown in figure (8-3A) 
and (8-3B) respectively. 

The exact transient and other characteristic 
values which were read from the computer print output have 
as follows: 
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(1) Step Input. 



Steady state output value 

Time to peak of first overshoot 

Settling time 

Percent overshoot 

Transient oscillatory frequency 



tp = 3.2 sec 
tg = 5.4 sec 
Mp = 10.2% 

= 0.92 rad/sec 



(2) Ramp Input . 



Steady state output value 
Steady state output rate 
Settling time 




R-0.002 
0.1 = R 
5.4 sec 



d. Fourth Modified Structure 



The fourth modified structure which is shown in 

figure (8-4) requires a feedback filter G (s) . After the 

s c 

* 

transfer function G sc (s) was evaluated from equations (8-27) 
and (8-36) the system was simulated in the computer. The 
transient response characteristics for a step R=1 Volt 
and a ramp R=0.1 rad/sec input, which were obtained in 
this case were the same with these corresponding to the 
first modified structure. They are shown in figure (8-4A) 
and (8-4B) respectively. 
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Fifth Modified Structure 



e . 



This modified structure involves a different 

If 

feedback singular filter G sc (s) whose transfer function 
was evaluated from equations (8-33) and (8-36). The 
transient response characteristics for a step R=1 Volt 
and a ramp R=0.1 rad/sec input, which were obtained 
when this structure was simulated in the computer are 
shown in figure (8-5A) and (8-5B) respectively, which are 
the same with these corresponding to the second modified 
structure . 



f. Sixth Modified Structure 

This is shown in figure (8-6) and was implemented 
from the previous one by introducing a gain adjustment as 
it is indicated in figure. 

Simulation of this system in the computer for 
a step R=1 Volt and a ramp R=0.1 rad/sec input, gave the 
same transient and steady state characteristics as in the 
case corresponding to the third modified structure. They 
are also shown in figures (8-6A) and (8-6B) respectively 
for each input which was considered. 
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5. Comment 



The above presented results which were obtained 
from the computer not only confirm the corresponding 
analytical results but also make clear the flexibility 
which the designer has when compensation is performed by 
using singular line theory. 
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IX. 



DESIGN OF COMPENSATION BY USING A 



SINGULAR CASCADE COMPENSATOR 



A. INTRODUCTION 



Consider the plant which is shown in figure (4-1), 
and that this plant is going to be compensated by a sin- 
gular cascade compensator. The open loop transfer function 
of this plant is: 



G (s) 



kFg (s) 



C s + z 1 ) (s+z 2 ) . . . (s + z m ) 
s N (s+Pj) (s+p 2 ) . . . (s+p n ) 



(9-1) 



It was assumed that the above transfer function is 
known and that the problem's specifications have been 
interpreted in terms of the desired dominant characteristic 
root (s, co n ) for the compensated singular system. It was 
also assumed that tolerance on the polar coordinates z, and 
oi n values for this root have been established from the 
problem's specifications. 

The compensated singular system is shown in figure 
(6-6) where G sc (s) denotes the transfer function of the 
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required singular cascade compensator, whose standard form 
is : 

(9-2) 



3 (s + 2 ) (s+Z ) . . . (s+Z ) 

G Cs) - ^ 

SC N 

S (s+Pj) (s+p 2 ) . . . (s+p n ) +a(s+z i ) (s+z 2 . . . (s+z m ) 

as was derived in Section VI C. 

Therefore the design problem in general will be limited 
to the problem of determining the required parameter pair 
values. For different types of applications this problem 
can be stated in terms of the parameter plane as: What must 
be the operating point, or the operating line segment, or 
the operating region of the parameter plane in order that 
the performance of the compensated singular system be as 
desired? 

The required design procedure in all these cases has 
been developed in Section VI. In this section the design 
process in the case where an operating point is desired to 
be defined on the parameter plane and from this the pre- 
cise form of the transfer function (9-2), was formulated 
in a sequence of basic design steps which must be followed 
in a case. 

B. DESIGN PROCESS 

The required design process in general when an ac- 
ceptable operating point of the parameter plane is desired 
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to be found, is relatively simple and can usually be 
completed with satisfactory accuracy for most applications 
with graphical and mathematical techniques, by the fol- 
lowing four basic design steps, 

1 . First Step 

The first step consists of drawing the Root locus 
of equation: 

kF (s) = -1 (9-3) 

& 

2 . Second Step 

The second step consists of insertion of an ini- 
tial cascade compensator, when the shape of the above Root 
locus is such that either an acceptable operating singular 
point cannot be located on this Root locus or although it 
can be located it is not acceptable under dominancy consid- 
erations (See VI, F,2). 

When an initial cascade compensator is required 
then this requirement will actually affect only the 
transfer function of the singular cascade compensator and 
neither the complexity of the rest of the design process 
nor the compensation structure (See example (6-6) and 
figure (6-28)). 

In the simple case where such an initial compen- 
sator will not be required then this step will consist of 
determining an acceptable operating singular point on the 
Root locus of equation (9-3), 
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3 . Third Step 



The third step consists of determining the 
equation of the singular line corresponding to the 
selected operating singular point (See VI D3a equation 
(6-29)). 

4 . Fourth Step 

This step which is the last one consists of de- 
termining the required operating point of the parameter 
plane, i.e., the required pair of a and g parameter values 
(See examples (6-6), (6-7), (6-8)). 
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X. 



EXTENSION OF THE SINGULAR LINE THEORY 



TO PARAMETER SPACE 



A. INTRODUCTION 



The parameter space is an extension of the parameter 
plane when there are three or more parameters. 

Considering the characteristic equation: 



n 

FCs) = H a, s k = 0 (10-1) 

k=0 K 



where n is a real integer and a^ is a real valued coef- 
ficient which may be a function of three or more parameters, 
i . e : 

a k = g( a ,3,Y,$, *•* nth ) (10-2) 

-L -L 

where n c denotes the n c parameter. 

The point s can be defined in polar coordinates as: 



-C 



am 



+ J 



to 



n 



/7 



(10*3) 



B. SINGULAR CONDITIONS 



If for a given pair u ns values, where [? S |<1 and 

a> ns >0, equation (10-1) has an infinite number of real valued 
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1" "h 

sets of a,e,y,<5, ... n parameter values which satisfy 
it, then for this equation singular conditions exist. In 
this case the point w ) is a singular point. This 

equation also for which singular conditions exist is 
called singular. 

C. SINGULAR CONDITIONS ANALYSIS IN THE LINEAR CASE 



A simple case which was analyzed is the linear case, 
in which the variable parameters appear in the coefficients 
of the characteristic polynomial in the general form: 

a k = b k a + c k 6 + d k Y + e k S + + ••• 

t h t h 

+{ Coefficient of the n parameters }^{n r parameter}+ 
+{Constant term}^ 



where b^, c^., ... are real constants. 

In this case equation (10-1) can be rearranged in the 



form: 

f 1 (s)a + f 2 (s) e + f 3 (s) y + . . . +f (s) {n tb parameter} + 

+ f n+1 Cs) - 0 (10-4) 

where f-^(s), f 2 (s), ... are polynomials of s. By analogy 
as in the case of the parameter plane (See Section IIB2) 
equation (10-4) can be written as: 



r i e ' 



j 0 i 



«+r 2 e' 



J 0 : 



6+r 3 e^ 9 3 



y+. 



‘ +r n e 



^" 9 n{n tb parameter }+r n+1 e^ 0n+1 



=0 
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or 

TjCcose 1 +jsin0 1 )a+r 2 Ccos9 2 +jsine 2 )B+. . .+ ... 

th 

+r n Ccos9 n +j sin9 n ) {n parameter }+r n+1 (cos9 n+1 +jsin9 n+1 ) = 0 

Equating then real and imaginary parts to zero, 
yields : 

r.sin9 ct+r sin9„ 3+... + r sin 9„{ n } + r ^ sin9 = 0 
1 1 2 2 n n n+i n+i 

and 

1"h 

r 1 cos9. ot+r cos9. B+...+ r cos9 { n } +r' cos9 = 0 
* i 2 2 n n n+in+i 

According to [1] the above system has an infinite 
number of solutions if and only if: 

r t sing r 2 sine 2 

r,cos0, r„cos9 

1 1 2 2 

or 

tanfl, = tan9 = tan9„ =...= tan9 = tan9 (10-5) 

l 2 3 n n+ 1 v J 

Equation (10-5) implies that for the characteristic 
equation (10-1) or C10~4) singular conditions exist if 
and only if: 



r 0 sin9, r sin9 r sine , 

3 3_ _ n n _ n+i n+i 

r cose r cos9„ r „x,cos9^. , 

3 3 n n n + 1 n+ i 
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/V s ) = /W s ) + K i" 

/V s ) ■ * K 2 " 

( 10 - 6 ) 

/V s ) = Ari (s) ' K J ' 



/V s ) "/W s )* K n n 

Assuming that f ^(s) ^ 0 and dividing equation 
(10-4) by £ n+1 (s) yields: 



^(s) 

Vi Cs) 



£,Cs) f,(s) 

a+ 6+ — 



f nV s ) 



w s ) 



V s ) 

W s > 



n th +l=0 



(10-7) 

or 



F 1 (s) a+ F 2 (s) e+ F 3 (s) y+...+ F n (s) n th +l = 0 

(10-8) 

When the characteristic equation is considered in the 
form of equation (10-8) then the necessary and sufficient 
conditions for singular conditions given by equation 
(10-6) become: 
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(.10-9) 



/F 1 (s ) = K x n 

/V s > = K z" 

Al£ S) * V 
An^'-'W 

where any of the K^, K 2> ..,K n is real integers. 

D. SINGULAR THEORY AND CONTROL SYSTEMS 

From equations (10-8) and (10-9) and by the same pro- 
cedure which was used in Section III for the case of the 
two parameters, it was found that in this case of the 
parameter space where there are n parameters in the 
characteristic equation of a control system the general 
form of a singular characteristic equation is: 

ni-t m ? ni . ■* m 

aC 1 F ± (s)+ 6C 2 F L (s) + y C 3 F 5 (s)+...+ {n tn }C n F n (s)+C n+1 =0 

(10-10) 

where C. , C~, C,,..., 0,0^, are real constants and m, , 
m 2 , m^,..., m n are n different real integers. (Otherwise 
when there are two or more of these identical then the 
number of the actual parameters can be reduced by combining 
terms of equal powers of F(s)). 
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For equation CIO - 10 ) the necessary and sufficient 
conditions, for singular conditions are given by equation: 



where K is a real integer. 

Equation (10-11) implies that all the points of the 
Root locus of equation: 



which do not lie on the real axis are singular points for 
the characteristic equation (10-10), where K in this case 
is considered to be the Root locus variable parameter which 
takes values from zero to plus infinity and from zero to 
minus infinity. 

Finally any control system with n variable parameters 
appearing in its characteristic equation linearly depen- 
dent, is singular, if and only if, its characteristic 
equation can be written in the form of equation (10-10) 
which in the simplest case takes the form: 



1 . Equation of a Singular n Dimensional Surface 

When an n dimensional parameter space was con- 
sidered then the general form of the corresponding singular 
characteristic equation of a control system is given by 
equation (10-10). Then according to the previous analysis 




( 10 - 11 ) 



KF(s) = -1 



( 10 - 12 ) 



1+ ctF(s) + BF 2 (s) + yF 3 (s) +...n th F n (s) = 0 (10-13) 
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when the value of a singular point s g is substituted in 
equation (10-10) it becomes: 

- - th ' 

aC^ + ( 5 C 2 + yCj +...+ n C n + = 0 (10-14) 

+ " >• + -* 

where C^, C 2 , C^,..., C n , C n+ ^ are real constants. 

Equation (10-14) is the equation of the singular n 
dimensional surface, corresponding to the singular point 
s s , which was considered. 

In order to determine a singular point of this 
singular surface the values of n-1 parameters must be 
defined and then the value of the parameter left can be 
evaluated from equation (10-14). Thus one among an infinite 
number of real valued sets of parameter values can be de- 
fined, which in other words defines the coordinates of a 
singular point in the parameter space. 

E. COMMENT 

From the previous analysis it can be seen that any 
non-linear case, either in the parameter plane on in the 
parameter space can be studied by defining any different 
parameter's product as a new variable parameter, e.g., the 
non-linear case in the parameter plane can be considered 
as a linear case in the parameter space (three dimensions 
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for this case) by defining the product ag equals with y 
which will be the third variable parameter. 

After transforming a non-linear characteristic 
equation in a linear characteristic equation as it was 
described above then all these which have been stated for 
the linear case can be applied for the transformed 
characteristic equation. 
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XI. 



CONCLUSIONS AND RECOMMENDATIONS 



FOR FURTHER WORK 



This thesis has developed practical methods for de- 
signing singular compensators, i.e., compensators which 
produce the desired singular conditions in a single input- 
single output control system. 

The key to the development of these methods was the 
form of the characteristic equation which was considered, 
i.e: 

af 1 (s)+Bf 2 (s)+f 3 (.s) = 0 

where: f-^(s), f 2 C s ) and f^Cs) are polynomials of s and 
a , B are real variable parameters. 

Necessary and sufficient conditions for the above 
characteristic equation to be singular were determined and 
the singular line theory was then investigated to make 
possible the- development of the required standard set of 
rules for the design engineer. 

Application of singular line theory in control system 
design led to the conclusion that it can be used in 
general for design of compensation of any single input 
single output linear control system with the following 
basic advantages. 
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Required filters for any desired compensation 
structure have transfer functions which are known as a 
function of the two parameters a and 3. The parameter 
plane diagram required for analysis and synthesis of a 
design problem consists of straight lines, i.e., the 
singular lines which can be drawn by hand, without needing 
a computer. Design techniques and procedures are charac- 
terized by simplicity, flexibility and accuracy. Finally, 
tolerances on the parameter a and g values can be pre- 
cisely established. 

When considering the sensitivity of the transient 
response to parameter variations, it was found that the 
control system is less sensitive to variations in a and g 
if the operating point on the singular line is remote from 
the point at which the singular line is tangent to the 
constant and w ns curves. 

The basic design tools of the singular line theory 
design method are: The Root locus of the plant and the 

parameter plane diagram of the compensated singular 
system, i.e., under certain conditions or in other words 
singular conditions, a classical design method (The Root 
locus method), and a modern design method (the parameter 
plane method) can work together in such a way that each 
compensates disadvantages of the other. 
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For example stability analysis which is not a simple 
process in the parameter plane is performed on the Root 
locus diagram in the s-plane or studying the effect of 
variations in either one or both of the parameter a and g 
values on the transient response of the system, which is 
not a simple process on the Root locus diagram is per- 
formed in the parameter plane, etc. 

Conceptually the compatibility of the Root locus 
classical design method, which is adequate when there is 
a single variable parameter with the parameter plane 
modern design method, which is adequate when there are 
two variable parameters, when singular conditions exist, 
is based on the fact that in this case one of the para- 
meters can be defined through the equation of a singular 
line, as a function of the other. Thus a two parameter 
problem becomes a one parameter problem. The similarity 
also of the Root locus of the plant with the Root locus 
of the compensated singular system, independently of which 
singular line is considered generalized the validity of 
the above concept and made possible the development of 
graphical design techniques which simplify the work and 
markedly increase the power of visualization of the de- 
signer . 

It was also realized that the singular line theory 
could be applicable in self-adaptive control system, since 
movement of the operating point (a,g) along a dominant 
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line segment of a singular line or within a dominant root 
region in the parameter plane implies almost fixed system 
performance. In such a case the adaptive controller need 
only to drive the system back to the dominant root line 
segment or dominant root region and not to a specific 
operating point in the parameter plane. 

The singular line theory was further extended in the 
parameter space and in the most general case where products 
of the variable parameters appear in the characteristic 
equation. 

Further study is needed for the case of a singular 
system with a negative error constant. This study can be 
initialized from an analysis of equation (8-25) . 

Also the correlation, if any, of a singular feedback 
compensator with the linear state variable feedback maybe 
interesting . 

Some other related areas for future work, which are 
recommended are: Applications of the singular line theory 
in self-adaptive control systems. Singular surface theory 
and control systems, singular line theory and multi-input, 
multi-output linear control systems, and finally any of 
the above in the general case where products of the variable 
parameters also appear in the characteristic equation, i.e., 
non-linear cases. 
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TABLE OF SINGULAR POINTS FOR EXAMPLE (5-1) 
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APPENDIX D 

TABLE OF ROOTS OF THE CHARACTERISTIC EQUATION FOR EXAMPLE (6-2) 

CORRESPONDING TO THE SINGULAR POINT ( **> = 0.40, a> n <rO. 97096) FOR DIFFERENT PARAMETER PAIRS 
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TABLE I 

FIRST SINGULAR CASE 
RESULTS OF OUTPUT AND ERROR ANALYSIS 
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TABLE II 

SECOND SINGULAR CASE 
RESULTS OF OUTPUT AND ERROR ANALYSIS 
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TABLE I 

EXAMPLE (6-5): SINGULAR POINT (C s=0. 40, u> ns =0. 97096) 

Transient Response Characteristic Values 
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System Unstable 



TABLE II 

EXAMPLE ( 6 - 5 ): SINGULAR POINT ( = 0 . 3 ,a) ns = 1 . 15268 ) 

Transient Response Characteristic Values 
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System Unstable 



TABLE III 

EXAMPLE (6-5): SINGULAR POINT ( C s = 0.75, ^ = 0.060667) 
Transient Response Characteristic Values 
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APPENDIX G 
TABLE I 

EXAMPLE (6-7): OPERATING POINTS AND CORRESPONDING SINGULAR 

POINTS AND SINGULAR LINES. 
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APPENDIX H 



SINGULAR POINT PROGRAM 



The purpose of this program is to determine the roots of a 
characteristic equation which are singular points with respect to the 
parameters a and $. 

This program is applicable to polynomials whose coefficients are 
of the form (b a+c&+d), where a and g are variable parameters and b, c, 
and d are real constants. 

For a given characteristic equation this program solves the 
system (2-4), by the Cramer's rule method for a desired value of t, 
and for a range of values for an, which has been determined within 
the program and can be changed when it is desired, by changing the 
card 0035. 

The four columns which are generated at the print output are the 
following: The first is a list of the values for a n , which were used 

for the calculations and the other three are the corresponding nume- 
rators and common denominator from the Cramer's rule solution of the 
system (2-4). 

Any value of «) n in the first column, for which the associated 
values of the other three columns are zero's, together with the value 
of t which was used defines a singular point for the characteristic 
equation. 

The input data for the program are: 

a. First card: 

Value of 4 and N, format (F4.2, 12) where C is the desired 
value of damping ratio and N is the order of the characteristic equation. 

b. Second card: 

Constant coefficients in ascending order (8E 10.5) 

c. Third card: 

Alfa coefficients in ascending order (8E 10.5) 

d. Fourth card: 

Beta coefficients in ascending order (8E 10.5) 
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APPENDIX I 



FIGURES' 




Block diagram of a unity feedback system 




H(s) 



R E 






FIGURE (4-2) 

First compensation structure 
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FIGURE (4-3) 

Second compensation structure 




FIGURE (4-4) 

Third compensation structure 
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FIGURE (4-5) 

Fourth compensation structure. 




FIGURE (4-6) 

Fifth compensation structure. 
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X- SCALE -2.00E+0I UNITS INCH. 
Y* SCALE' 6 OOEK)l UNITS INCH. 



FIGURE (5-1) 

Example (5-1) :Parameter plane 



di agram. 
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\ 



X-SCALE =2.00 E+01 UNITS INCH. 

Y-SCALE =6.00 E+01 UNITS INCH. 

FIGURE (5-2) 

Example ( 5-1 ): Parameter plane diagram showing 
a singular line. 
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FIGURE (6-1) 

Block diagram corresponding to the second singular 




X-SCALE- 5.00E-OI UNITS INCH. 
Y- SC ALE =4 OOEKX) UNITS INCH. 



FIGURE (6-2) 

Example ( 6-1 ): Parameter plane diagram. 
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FIGURE (6-3) 

Block diagram of the plant. 




FIGURE (6-4) 

Block diagram of the compensated singu- 
lar system corresponding to the "second 
singular case" ,s cheme . 
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FIGURE (6-5) 

Block diagram of the compensated singular 
system corres pondi ng to the modified second 
singular case, scheme. 




FIGURE (6-6) 

Simplified block diagram of the modified 
second singular case, scheme. 
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FIGURE (6-7) 

Example (6-2): Root locus of the uncompen- 
sated system. 
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X-SCALE- I.00E+02 UNITS INCH. 

Y- SCALE 3 l.0 0E-*-02 UNITS INCH. 

FIGURE (6-8) 

Example (6-2): Parameter plane di agram,K=20 . 
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FIGURE (6-9) 

Example (6-2): Root locus of the compensated singu- 
lar system corresponding to the singular root: 
s = -0.3884±j 0.8899 or ( C $ =0 . 40 ,“ ns =0 . 97096 ) 
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X-SCALE » 1.00 E-*-02 UNITS INCH 
Y-SCALE » 1.00 Ef02 UNITS INCH. 



FIGURE (6-10) 

Example (6-2): Parameter plane di agram,K=35. 
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FIGURE (6-1 1 ) 

Example (6-3): S-plane diagram. 
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FIGURE (6-12) 

Example (6-4): S-plane diagram. 
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FIGURE (6-1 3a ) 

Example (6-4): Introduction of an additional zero. 
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FIGURE (6- 1 3b ) 

Example (6-4): Introduction of an additional pole. 
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FIGURE ( 6 - 1 3 c ) 

Example (6-4): Introduction of an additional pole and zero 
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FIGURE (6-1 4a ) 

Example (6-5): Block diagram of the compensated sin- 
gular system with the singular compensator considered 
i rr i nt eg rated form . 



s ( s + 1 ) ( s + 5 ) I G( s )“ s ( s+1 ) (s + 5 ) 
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FIGURE (6-1 4b ) 

Example (6-5): Block diagram of the compensated sin- 
gular system with the singular compensator considered 
to consist of distinct components. 
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V-Scal e=l . 00 Units/Inch. 

FIGURE (6-17) 

Example (6-5) :Time response of the compensated singular system for a step input 
(Parameter values: A=20 ,B=3 . 712) r r 
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FIGURE (6-19) 

Examp 1 e ( 6 - 5 ) : Root locus of the compensated singular 
system corresponding to the singular root (cc=0.30, 
(u ns =1 .1 5268 ) . 
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FIGURE (6-24) 

Example (6-5): Root locus of the compensated sin- 
gular system corresponding to the singular root 
U s =0 .75, a>n s = 0 .60667) . 
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FIGURE (6-26) 

Example (6-6): Block diagram of 
the plant. 
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FIGURE (6-27) 

Example (6-6): Block diagram of the 
initially compensated plant. 



243 



B(s + 2) r / _ \ _ s+2 






CJ 



CM 



3 

+ 

O 

CM 

+ 

to 



LO 

+ 

to 



II 

to 

CJ 

to 

CD 



L> 



oc 



E 

CD 

-M 

to 

>> 

to 

S- 

03 

3 

CD 

C 

•r- 

tO 

T3 

<D 

03 
to 
C 
CD 
CL 
03 E 
CO O 
CM O 
I 

tO O) 



a: 4- 

=> o 

co 

— E 

LL_ 03 

s- 

CD 

03 



CJ 

o 

CO 



to 

I 

to 



CD 

a. 



03 

X 



244 



e(s + 2) ^ s + Z 20 

s ( s + 1 ) (s + 5) (s+20f f s + 20 G l s ) “s ( s + 1 ) ( s + 5 ) 



^ , 




CD 

+-> 

i/) 

CO fO 

CO 
s- CM 
fO I 

I— <£) 



CO 

= c d 

•r- S_ 
CO 3 
CO 
"O *p— 
CD 4- 



-Q 

CO 

CM 

I 

CO 



fC c 

CO *r- 

c 

CD cj 

E o 

o jc 

u to 



CD +-> 
to 

+-> JC 
UJ +-> 

DC 4— 

=) Or 

CD +_> 



ll fo 5 

S- 

03 E 
fO n3 
*•“ S- 
"O CO 
CO 

-X -r- 

O -o 
O 

' _s^ 

CD O 
O 



^ -Q 
CO 
I +-> 

c 

— " CD 
r— 
CD (O 
1 — > 
CL.r- 

E =3 
co cr 

x UJ 
UJ 



245 




I 

d -d oo 

cnu 

d -r- O 
•r- _d +-> 
w 2 d 
uo 

“O s- c 

CD O CD 
+->+-> Q_ 
rd fd £ 
CO CO O 
d d o 
a) aj 

Cl a_r— 

£ E ro 

O O -r- 

o o +-> 



a) r— c 

_d r— ’I — 

+-> fa 
S- CD 
4— CL) _d 

0 > +-> 

^ o 

0 E 

CO to CD d 

oo i. -d to 

1 cn +-> 

CO H 3 £. 

- — " * i — co to 

“O +-> i — 
LjJ C 3 
CD CT 
D U W C 
CD O CD * i — 
*— < i — i- CO 
Ll_ _q q_ 

CD CD 
+-> S- _d 
d +-> 
CD * 

i — LD CD 
fO — - d 
> O- 
•r- COC 
d CD v- 

cr _q 
LD CD E 

i- o 
••(DU 
'-d 

1 _Q 

E 

— " CD “O 
+-> CD 

CD co E 

r- >> 

Q. l/l O 
E < 4 - 
(0 i- 
x <o oo 

LU i — •i—* 



246 




FIGURE (6-29) 

Example (6-6): Root locus of the uncompensated system. 
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FIGURE (6-30) 

Example (6-6): Root locus of the plant after 
it was compensated by an initial cascade com- 
pensated with transfer f unc ti on ( s + 2 )/ ( s + 20 ) . 



248 




FIGURE (6-31 ) 

Example (6-6):Root locus of the singular 
compensated system corresponding to the 
singular root ( r =0 .6 , u» =1 .2 326 ) or (s = 
-0 .7395±j0 .9861 ) s ns 
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r i / c c \ FIGURE (6-32) 

(Parameter'val uesf A^^°B = 28?987) e compensated s1 "9 u,ar ^ a step input. 
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FIGURE (6-37) 

Example (6-6):Time response of the compensated singular system for a step input. 
(Parameter values: A=131.5, B=216. 141 ) 
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Y-Scale=1.00 Units/Inch. 

FIGURE (6-46) 

Example (6-6):Time response of the compensated singular system for a step input. 
(Parameter values: A=1700 , B = 541 1 . 722 ) 
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FIGURE (6-47) 

Example (6-6):Time response of the compensated singular system for a step input. 
(Parameter values: A=1800, B=5742.968) 
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o/o 800 joo boo zoo \ ooo I zvo 




X-SCALE=5.00E+01 UNITS INCH. 

Y-SCALE = 2.00E + 02 UKlITS INCH. 

FIGURE (6-49) 

Example (6-6): Parameter plane diagram of the 
compensated singular system. Si ngul ar point 
U s =0. 6, w ns = l. 232586) 
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X-SCALE=4. 00E+02 UNITS INCH. 
Y-SCALE=1 . 00E+03 UNITS INCH. 

FIGURE (6-50) 

Example (6-6): Parameter plane diagram of the 
compensated singular sy s tern . S i ngu 1 ar point 
(? s = 0.6,(D ns = 1.2325) 



268 







FIGURE (6-51 ) 

Example (6-7): Root locus of the compensa- 
ted singular system. 
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FIGURE (6-52) 

Example(6-7):Para meter plane diagram and 
bo un da ry lines. 
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FIGURE (6-53) 

Example (6-7):Time response of the compensated singular system for a step input. 
(Parameter values: A=900, B=2400) 
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FIGURE (7-1) 

Example (7-1): Block diagram of the 
i ni ti al sys tern (Plant) . 




FIGURE (7-2) 

Example (7-2): Block diagram of the 
compensated singular system. 
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FIGURE (7-3) 

Root locus diagram. 
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FIGURE (7-4) 

Root locus of the compensa- 
ted singular system. 
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First modified singular compensation 
structure. 




structure. 
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FIGURE (8-3) 

Third modified singular compensation 
structure . 




FIGURE (8-4) 

Fourth modified singular compensation 
structure. 
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FIGURE (8-5) 

Fifth modified singular compensation 
structure. 




FIGURE (8-6) 

Sixth modified singular compensation 
structure . 
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-Scale=5.00 Units/ inch 
-Scale=0.40 Units/inch 
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FIGURE ( 8-3A) 

Third modified singular s tructure . Ti me response of the system for a step input. 
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FIGURE ( 8-4A) 

Fourth modified singular structure. Time response of the system for a step input. 




ZJ 

Q_ 

C 



03 

s- 



03 



s- 

o 



O) 



CO 

>> 

CO 



O) 



<1) 

CO 

o 

^ Q_ 

QQ CO 
d) 

I s- 

00 

^ d) 

E 

LxJ -r- 
CC h- 



CD <1> 
i— i S- 
Ll_ 13 



U 

13 

S- 

+-> 

CO 

s- 



X >- LI- 



289 



-Scale=5.00 Units/ inch 
-Scale=2.00 Units/ inch 
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